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ABSTRACT 

We present measurements of galaxy clustering from the Baryon Oscillation Spectroscopic 
Survey (BOSS), which is part of the Sloan Digital Sky Survey III (SDSS-III). These use the 
Data Release 9 (DR9) CMASS sample, which contains 264 283 massive galaxies covering 
3275 square degrees with an effective redshift z = 0.57 and redshift range 0.43 < z < 
0.7. Assuming a concordance ACDM cosmological model, this sample covers an effective 
volume of 2.2 Gpc 3 , and represents the largest sample of the Universe ever surveyed at this 
density, n sw 3 X 10~ 4 /i _3 Mpc 3 . We measure the angle-averaged galaxy correlation function 
and power spectrum, including density-field reconstruction of the baryon acoustic oscillation 
(BAO) feature. The acoustic features are detected at a significance of 5<7 in both the correlation 
function and power spectrum. Combining with the SDSS-II Luminous Red Galaxy Sample, 
the detection significance increases to 6.7a. Fitting for the position of the acoustic features 
measures the distance to z — 0.57 relative to the sound horizon Dy jr s = 13.67±0.22 at z = 
0.57. Assuming a fiducial sound horizon of 153.19Mpc, which matches cosmic microwave 
background constraints, this corresponds to a distance Dy(z — 0.57) = 2094 ± 34Mpc. At 
1 .7 per cent, this is the most precise distance constraint ever obtained from a galaxy survey. 
We place this result alongside previous BAO measurements in a cosmological distance ladder 
and find excellent agreement with the current supernova measurements. We use these distance 
measurements to constrain various cosmological models, finding continuing support for a flat 
Universe with a cosmological constant. 

Key words: cosmology: observations, distance scale, large-scale structure 



1 INTRODUCTION 

Explaining the late-time acceleration of the expansion rate of the 
Universe (Riess et ali. 1998; Perlmutter et al. 1999) is one of the 
most perplexing problems in modern physics. All known attempts 
require exotic ingredients: a new, very small energy scale in a cos- 
mological constant or low-mass field, a change to General Rela- 
tivity to weaken gravity on large scales or at low densities, or ex- 
tra dimensions of space-time. Empirical observations will provide 
clues as to the cause by providing precision measurements of the 
expansion history and the growth of cosmological structure over 
time (e.g. |Kolbetal.|2006> . 

One of the key methods for measuring the expansion history 
is to use features in the clustering of galaxies within galaxy surveys 
as a ruler with which to measure the distance-redshift relation. Ob- 
taining precision distance measurements is a long-standing chal- 
lenge in astronomy, and the Baryon Acoustic Oscillations (BAO) 
signal in the two-point clustering of galaxies provides a particularly 
robust quantity to measure. The BAO arise because the coupling of 
baryons and photons by Thomson scattering in the early universe 
allows acoustic oscillations at early times, which in turn leads to 
a rich structure in the distribution of matter and the anisotropies 
of the cosmic microwave background (CMB) radiation. The dis- 
tance that acoustic waves can propagate in the first million years of 
the universe becomes a characteristic comoving scale (Peebles & 



|Yu|1970[|Sunyaev & Zel'dovich|1970l|Doroshkevich et al.|1978l 

a description of the physics leading to the features can be found in 
Eisenstein & Hu 1998 or Appendix A of Meiksin, White & Peacock 
1999 and a discussion of the acoustic signal in configuration space 
can be found in Eisenstein et al. 2007b). As the acoustic signature 
is imprinted on very large scales (~ 150 Mpc) it is quite insensi- 
tive to astrophysical processing that occurs on smaller scales, thus 
BAO experiments are affected by a very low level of systematics 
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induced by such processes. Recent reviews of BAO as a probe of 
dark energy may be found in Eisenstein ( 2008 1 and Weinberg et al. 
l |20i"2l 

This acoustic signature has now been detected in many dif- 
ferent galaxy surveys, using a variety of methods to analyse the 
evolved density field (Percival et al. 2001; Miller, Nichol & B 
tuskil[200Tl pisenstein et al.||2005UCole et al.||2005HHutsi||200r| 
Padmanabhan et al. 2007, Blake et al. 2007, Per cival et al.|2 007c; 
Okumura et al. 2008 , Gaztanaga, Cabre & Hui 2009 ; Kazin et 
2010; Percival et al. 2010; Reid et al. 2010, Blake et al. 201 



Beutler et al.|201 l[|Seo et al.|2012 , Padmanabhan et al. 2012a}, and 
it is already producing stringent constraints on cosmological mod- 
els. Constraints can be obtained from either photometric or spectro- 
scopic samples, though for the same volume and number of galax- 
ies the spectroscopic samples provide much stronger constraints. 
The first BAO measurements came from the 2dF Galaxy Redshift 
Survey (2dFG RS; |Colless et al|2003) and the Sloan Digital Sky 
Survey (SDSS; York et al. 2000); when combined, the most recent 
analyses give a 2.7 per cent measurement of the distance-redshift 
relation at z = 0.275 (e.g. Percival et al. 2010). Adding to these 
data, |Blake et al.| J201 la| ) measured the BAO feature at z = 0.6, 
using the WiggleZ survey (Drinkwater et al 2010), making a 4 per 
cent distance measurement from 132 509 galaxies. This result was 
subsequently improved to provide distance measurements of accu- 
racy 7.2 per cent, 4.5 per cent and 5.0 per cent in three bins centered 
at redshifts z = 0.44, 0.60, and 0.73 respectively, using the full 
sample of 158 741 galaxies from this survey {Blake et al.|2 01 lb). 
|Beutler et al.| ( |20l"T| ) made a 4.5 per cent measurement at z = 0.1 
with the 6dF Galaxy Redshift Survey (6dFGRS: [Tones et aL|2009) . 
Thus the BAO technique has recently provided a distance-redshift 
relation at a series of redshifts both higher and lower than the 2dF- 
GRS and SDSS measurements. 

The BAO signal in an evolved galaxy field, such as those anal- 
ysed in the papers described above, differs from that predicted in 
the matter field by linear theory alone. The dominant difference 
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is caused by matter flows and peculiar velocities on intermediate 
scales (~ 20 h" 1 Mpc), which act to suppress small-scale oscilla- 
tions in the galaxy power spectrum and smooth the BAO feature in 
the correlation function (Eisenstein et al. 2007b, Crocce & Scoc- 
|cimarro|2008j |Matsubara|2008a|b| >. |Eisenstein et al.| ( |2007a) sug- 
gested that this smoothing can be reversed, in effect using the phase 
information within the density field to reconstruct linear behaviour. 
Although not a new idea (e.g. Peebles 1989 . 1990 , Nusser & Dekel 
1992; Gramann 19931, the dramatic effect on BAO recovery had 
not been previously realised, and the majority of the benefit was 
shown to be recovered from a simple reconstruction prescription. 
This reconstruction technique has been used to sharpen the BAO 
feature and improve distance constraints on mock data ( Padmanab- 
han & White|2009||Noh, White & Padmanabhan|20091 |Seo et"aT| 



2010||Mehtaetal.|2011^ , and it was recently applied to the SDSS-II 
Luminous Red Galaxy (LRG) sample (Padmanabhan e t al.|2012a} . 
The reconstruction was particularly effective in this case, provid- 
ing a 1.9 per cent distance measurement at z — 0.35, decreasing 
the error by a factor of 1.7 compared with the pre-reconstruction 
measurement. 

This study is the first in a set of papers to describe the clus- 
tering of galaxies at z ~ 0.6 from Data Release 9 (DR9) of 
the Baryon Oscillation Spectroscopic Survey (BOSS; Schlegel 



White & Eisenstein 2009), which is part of SDSS-III (Eisenstein 
et al. 2011). We present cosmological results based on fits to the 
BAO signature in the clustering of 264 283 galaxies in this paper. 



Redshift-Space Distortion (RSD) and Alcock-Paczynski lAlcock 
& Paczynski|1979| AP) measurements are presented in Reid et al. 



(2012JI, and an interpretation of these results in terms of dark en- 



ergy and modified gravity models is presented in Samushia et al. 
(2012). |Tojeiro et aL] ([20l2a) describe a new method for improving 
RSD measurements. Further constraints from fitting models to the 
full shape of the correlation function are presented in |Sanchezetal.| 
(2012). [Nuza et af] (201 21 have compared the clustering with the 
outcome of a large-volume cosmological simulation. 

Each of these papers focusses on the high redshift galaxy 
sample from BOSS, denoted "CMASS", where a set of colour- 
magnitude cuts are used to select a roughly volume-limited sam- 
ple of massive, luminous galaxies from a redshift of 0.43 to 0.7. 
We describe the construction of the galaxy catalog and measure- 
ment of both the correlation function and power spectrum of this 
sample, before and after applying the reconstruction algorithm of 
Eisenstein et al. (2007a). We present the results of two pipelines 
for the analysis of BAO: one utilising the correlation function, and 
one utilising the power spectrum. While both statistics contain the 
same information, they often have different advantages in the usage 
and application in the literature. We compare and contrast measure- 
ments made on mock catalogues using both techniques, and apply 
both to measure and analyse the BAO distance scale using CMASS 
data. In companion papers we present weights used to correct for 
artificial density fluctuations caused by observational effects (Ross 
|et al.|2012 1 and a set of mock catalogues used to estimate statistical 
errors (Manera et al. 2012). An analysis of a lower redshift sample 
of galaxies from BOSS will be presented in Parejko et a!.| ( [20T2) >. 
Clustering measurements from a smaller sample of CMASS galax- 
ies (the first six months of BOSS data) were presented by |White"et| 
[aTT| l [20 1 1) > and used to constrain halo occupation distributions, but 
these measurements did not extend to the BAO scale. 

The layout of our paper is as follows. We introduce the data in 
Section[2]and the catalogue used in Section[3] Analysis techniques 
are described in SectionH] and correlation function and power spec- 
trum measurements are described and presented in Sections Bland 



[6] respectively. These are compared and our final distance mea- 
surement presented in Section [7] This measurement is placed in a 
cosmological context in Sections [8] and [9] A brief discussion is 
given in SectionflO] Finally, a series of Appendices test the validity 
of various aspects of the methods used. 

Throughout the paper we assume a fiducial ACDM+GRHflat 
cosmological model with Q m = 0.274, h = 0.7, tt b h 2 = 0.0224, 
n s - 0.95 and a 8 = 0.8, similar to the best-fit WMAP 7-year 
model (Komatsu et al. 201 1 ). These parameters allow us to translate 
angles and redshifts into distances and provide a reference against 
which we measure distances. The BAO measurement allows us to 
constrain changes in the distance scale relative to that predicted by 
this fiducial model. 



2 THE DATA 



The Sloan Digital Sky Survey (SDSS; |York et al.pOOO) mapped 
over one quarter of the sky using the dedicated 2.5-m Sloan Tele- 
scope (Gunn et a l.|2006| ) located at Apache Point Observatory in 
New Mexico. A drift-scanning mosaic CCD camera (Gunn et al 



1998) imaged the sky in five photometric bandpasses (Fukugita et 
al.|1996||Smith et al.|2002||Doi et al.|2010| > to a limiting magnitude 
of r ~ 22.5. The imaging data were processed through a series 
of pipelines that perform astrometric calibration {Pier et aI7[2 003 ), 
photometric reduction (Lupton et al. 2001), and photometric cali- 
bration (Padmanabhan et al. 2008i. The magnitudes were corrected 
for Galactic extinction using the maps of Schlegel, Finkbeiner & 
Davis ( 19981 BOSS, as part of the SDSS-III survey ( [Eisenstein et 
al. 201 1 1, has imaged an additional 3100 square degrees of sky over 



that of SDSS-II (Abazajian et al. 2009) in the South Galactic sky, 
in a manner identical to the original SDSS imaging. This increased 
the total imaging SDSS footprint to 14 055 square degrees, with 
7600 square degrees at |fo| > 20 deg in the North Galactic Cap and 
3000 square degrees at |6| > 20 deg in the South Galactic Cap. All 
of the imaging was re-processed as part of SDSS Data Release 8 
( |Aiharaetal.|2011) . 

BOSS is primarily a spectroscopic survey, which is designed 
to obtain spectra and redshifts for 1.35 million galaxies over an ex- 
tragalactic footprint covering 10 000 square degrees. These galax- 
ies are selected from the SDSS imaging and are being observed 
together with 160 000 quasars and approximately 100 000 ancil- 
lary targets. The targets are assigned to tiles of diameter 3° using 
a tiling algorithm that is adaptive to the density of targets on the 
sky ([Blanton et al. 2003). Aluminium plates are drilled with 1000 
holes whose positions correspond to the positions of objects on 
each tile, which are manually plugged with optical fibres that feed 
a pair of double spectrographs. The double-armed BOSS spectro- 
graphs are significantly upgraded from those used by SDSS-I/II, 
covering the wavelength range 3600 A to 10 000 A with a resolv- 
ing power of 1500 to 2600 (Gunn et al. 2012). In addition to ex- 
panding the wavelength coverage from the SDSS-I range of 3850 to 
9220 A, the throughputs have been increased with new CCDs, grat- 
ings, and improved optical elements, and the 640-fibre cartridges 
with 3" apertures have been replaced with 1000-fibre cartridges 
with 2" apertures. Each observation is performed in a series of 
900-second exposures, integrating until a minimum signal-to-noise 
ratio is achieved for the faint galaxy targets. This ensures a homo- 
geneous data set with a high redshift completeness of > 97 per 



1 We do not consider any modifications to general relativity in this paper. 
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cent over the full survey footprint. A summary of the survey design 
appears in Eisenstein et al. (2011), and a full description will be 
provided in Dawson et al. (2012). 

2.1 Galaxy target selection 

BOSS makes use of luminous galaxies selected from the multi- 
colour SDSS imaging to probe large-scale structure at intermedi- 
ate redshift (0.2 < z < 0.7). The target selection is an extension 
of the targeting algorithms for the SDSS-II (Eisenstein et al.|2001) 
and 2SLAQ (Cannon et al. 2006) Luminous Red Galaxies (LRGs), 
targeting fainter and bluer galaxies in order to achieve the a number 
density of 3 x 10~ 4 /i 3 Mpc~ 3 . The majority of the galaxies have 
old stellar systems whose prominent 4000 A break makes them rel- 
atively easy to target using multi-colour data. The details of the 
target selection algorithm will be presented in Padmanabhan et al. 
(2012b); we summarise the details relevant to this paper below. 

The galaxy target selection in BOSS is divided into two 
classes of galaxies: LOWZ galaxies (0.2 < z < 0.43) and CMASS 
galaxies (0.43 < z < 0.7), analogous to the Cut-I and II SDSS-II 
LRGs. The 4000 A break resides primarily in the g and r bands for 
the LOWZ and CMASS redshift ranges respectively. The LOWZ 
sample in DR9 was somewhat compromised by a target selection 
error, now fixed, in the early data, and regardless it would have 
fewer objects and a lower effective volume than the SDSS-II LRGs 
over the same redshift range. We therefore restrict our analysis here 
to the CMASS sample and use the results from Padmanabhan et al. 
(2012a) for measurements in the lower redshift range. The small 
scale clustering results of the LOWZ sample is described in the 
companion paper of Parej ko et al.| ( |20i2) . Future BOSS analyses 
will use both the LOWZ and CMASS samples. 

The CMASS sample was designed to loosely follow a constant 
stellar mass cut (hence the name ConstantMASS) based on the pas- 
sive galaxy template of Maraston et al. (2009), and was designed 
to produce a uniform mass distribution at all redshifts. The distri- 
bution of CMASS stellar masses (Maraston et al. 2012) and veloc- 
ity dispersions (Thom as" et al.|2012| ) in various redshift bins con- 
firm that this goal was achieved. Unlike SDSS-II LRGs, we do not 
exclusively target intrinsically red galaxies with the CMASS cut. 
In fact, |Masters et al.| ( |20TT] showed that 26 per cent of CMASS 
galaxies are massive spirals. Therefore, whereas both the LOWZ 
and CMASS samples are colour-selected, CMASS is a significantly 
more complete sample than LOWZ at high stellar masses. This is- 
sue is discussed in detail in Tojeiro et al. (2012a), which consid- 
ers the passive evolution of galaxies between the SDSS-II Lumi- 
nous Red Galaxies (which form a subset of the LOWZ sample) and 
the CMASS sample. Most CMASS objects are central galaxies re- 
siding in dark matter halos of 10 13 h _1 M , but a non-negligible 
fraction are satellites that live primarily in halos about 10 times 
more massive (White et al. 201 1 Nuza et al. 2012). Galaxies in the 
CMASS sample are highly biased (b ~ 2), and bright enough to be 
used to trace a large cosmological volume with sufficient number 
density to ensure that shot-noise is not a dominant contributor to the 
statistical error in BAO measurements. The combination of large 
volume, high bias, and reasonable space density makes CMASS 
galaxies particularly powerful for probing statistical properties of 
large-scale structure. 

The CMASS target selection makes use of four definitions 
of flux computed by the photometric pipeline. All magnitudes 
have been photometrically calibrated using the uber-calibration of 
Padmanabh an et aT| (|2008) and corrected for Galactic extinction 
(Schlegel, Finkbeiner & Davis 1998). The model fluxes are com- 



puted using either a PSF-convolved exponential or de Vaucouleurs 
light profile fit to the r-band only, and are denoted with the mod 
subscript. Cmodel fluxes are computed using the best-fit linear 
combination of an exponential and a de Vaucouleurs light profile 
fit to each photometric band independently ( Abazajian et al. 2009), 
and are denoted with the subscript cmod. Point-spread function 
(PSF) fluxes are computed by fitting a PSF model to the galaxy, 
and are denoted with the subscript psf (Stoughton et al. 2002). Fi- 
bre fluxes are computed within a 2 arcsec aperture after the image 
is convolved with a kernel to produce a 2 arcsec FWHM PSF, and 
are denoted with the subscript fib2. Colours are computed using 
model fluxes. Magnitude cuts are performed on cmodel and fibre 
fluxes. 

The CMASS algorithm selects luminous galaxies at z > 0.4, 
extending Cut II of Eisenstein et al. (2001 1 to both fainter and bluer 
galaxies. We first select objects classified as galaxies by the imag- 
ing pipeline. These must then pass the following criteria: 



17.5 < i cmod < 19.9, 

'mod 'mod \ ^.U, 

d ± > 0.55, 

iflb2 < 21.5, 

icmod < 19.86 + 1.6(da. - 0.; 



(1) 

(2) 
(3) 
(4) 
(5) 



where the auxiliary colour d± is defined as {Cannon et al.| 2006) 

d± — r mod - i mod - (p mod - r mod )/8.0. (6) 

CMASS objects must also pass the following star-galaxy sep- 
aration cuts 

Vf-W > 0.2 + 0.2(20.0 -w), (7) 

z ps f - Zmod > 9.125 - 0.46z mod , (8) 

unless they also pass the LOWZ criteria (see Ross et al. 2012, Pare- 
[jko et al.|2012) , which only uses the default SDSS-II star-galaxy 
separation criterion. Stars will have essentially identical model and 
PSF fluxes, and this star-galaxy separation cuts on the difference 
between these two magnitudes. The slope with apparent magnitude 
in equations [7] and [8] which is not used in the standard star-galaxy 
separator of the photometric pipeline ( Strau ss et al.|2002] >, was set 
empirically by analysing commissioning spectroscopic data that re- 
laxed these cuts. Our choices yield a sample with approximately 3 
per cent stellar contamination, and it discards approximately 1 per 
cent of genuine galaxy targets, mostly at the faint end. The star- 
galaxy separation is known to fail when the seeing is poor, as PSF 
and model magnitudes approach one another for all object types in 
poor seeing. However this has been shown to have negligible effect 
on the angular distribution of targets in SDSS (Ross et al.|2011[ ). 

2.2 Masks 



We use the MANGLE software (Swan son et al.|2008 l to track the 
areas covered by the BOSS survey and the angular completeness 
of those regions. The mask is constructed of spherical polygons, 
which form the base unit for the geometrical decomposition of the 
sky. The angular mask of the survey is formed from the intersec- 
tion of the imaging boundaries (expressed as a set of polygons) and 
spectroscopic sectors (areas of the sky covered by a unique set of 
spectroscopic tiles) (see Blanton et al. 2003, Tegmark et al. 2004; 
|Aiharaetal.|2011| ). 

In addition to tracking the outline of the survey region and the 
position of the spectroscopic plates, we apply several "vetos" in 



© 201 1 RAS, MNRAS 000,pp3 



BAO in SDSS-III BOSS DR9 galaxies 




240° 220° 



200° 180° 160° 140° 
RA (degrees) 





^MH 






B 




0.7 0.8 0.9 I.C 


+20° 




completeness 


S> +10° --. 

"a 


« • 


i 


0° ■ 




m n y*" JB^ H 


-10° • 







0° 
RA (degrees) 



Figure 1. The sky coverage of the galaxies used in this analysis. The light 
grey shaded region shows the expected total footprint of the survey, totalling 
10 269 deg 2 . The coloured and dark grey regions indicate the DR9 spectro- 
scopic coverage of the survey, totalling 3792 deg 2 . Colours indicate the 
completeness within each sector used to build the random catalog as de- 
fined in Eq. |10| Sectors coloured dark grey are removed from the analysis 
by the cuts described in Section |3~5| The total effective area (accounting for 
all applied cuts and the completeness in every sector included) used in our 
analysis is 3275 deg 2 . The low completeness at many edges is due to un- 
observed tiles that will overlap the current geometry in future data releases. 



constructing the catalog. Regions were masked where the imaging 
was unphotometric, the PSF modelling failed, the imaging reduc- 
tion pipeline timed out (usually due to too many blended objects in 
a single field), or the image was identified as having critical prob- 
lems. Small regions around the centre posts of the plates where fi- 
bres cannot be placed due to physical limitations and around bright 
stars in the Tycho catalog (H0g et al. 2000) were also masked. The 
mask radius for stars from the Tycho catalog was 



R = (0.0802B 2 - 1.860B + 11.625) arcmin , 



(9) 



where B is the Tycho Bt magnitude clipped to fall in the range 
[6, 11.5]. We also place a mask at the locations of objects with 
higher priority (mostly high-z quasars) than galaxies. A galaxy can- 
not be observed at a location within the fibre collision radius of 
these points. In total we masked ~ 5 per cent of the area cover by 
the set of observed tiles due to our "veto" mask. The sky coverage 
of the CMASS galaxies is shown in Fig. [I] and the basic parameters 
including areas and galaxy numbers are presented in TablefT] 
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Figure 2. The galaxy number density as a function of redshift for the BOSS 
DR9 CMASS sample (thick blue line) used in this analysis, which ranges in 
redshift between 0.43 < z < 0.7. For comparison, we also plot the density 
for a SDSS-II DR7 LRG sample (thin red line) covering 0.16 < z < OAT, 
which was used in Padmanabhan et al. (2012al. Note that both selections 
include a small fraction of objects that fall outside the redshift cuts shown 
here. 



idlspec2ci software packager! The classification and redshift 
of each object are determined by fitting their coadded spectra to 
a set of galaxy, quasar, and star eigentemplates. The fit includes a 
polynomial background (quadratic for galaxies, quasars, and cata- 
clysmic variable stars; cubic for all other stars) to allow for residual 
extinction effects or broadband continua not otherwise described by 
the templates. The reduced \ 2 versus redshift is mapped in redshift 
steps corresponding to the logarithmic pixel scale of the spectra, 
Alog 10 (A) = 0.0001. Galaxy templates are fit from z — -0.01 
to 1.00, quasar templates are fit from z — 0.0033 to 7.00, and star 
templates are fit from z = -0.004 to 0.004 (±1200 kms -1 ). The 
template fit with the best reduced \ 2 is selected as the classification 
and redshift, with warning flags set for poor wavelength coverage, 
broken/dropped and sky-target fibres, and best fits which are within 
Ax 2 /dof = 0.01 of the next best fit (comparing only to fits with 
a velocity difference of less than 1000 kms -1 ). This method is the 
same as used for the SDSS DR8; see Aihara et al. ( 201 1 ) for further 
explanation. 

For galaxy targets, a dominant source of false identifications 
is due to quasar templates with unphysical fit parameters, e.g., large 
negative parameters causing a quasar template emission feature to 
fit a galaxy absorption feature. Thus, for galaxy targets, the best 
classification and redshift are selected only from the fits to galaxy 
and star templates, and we restrict to fits the pipeline classifies as 
robust. 

Fig. [2] shows the galaxy number density of the CMASS sam- 
ple, compared with the SDSS-II LRG sample. The CMASS galax- 
ies have approximately three times the density of the SDSS-II LRG 
sample, and sample the underlying density field with lower noise 
and higher fidelity. Although redshifts are recovered at higher and 
lower redshifts, we limit the CMASS redshift range to 0.43 < z < 



2.3 Measuring galaxy redshifts 

Spectroscopic calibration, extraction, classification, and red- 
shift analysis were carried out using the v5_4_45 tag of the 



http : //www . sdss3 . org/svn/repo/idlspec2d/tags/v5_4_4 5/ 
3 These fits are stored in the "*_NOQSO" versions of the Z, Z_ERR, 
ZWARNING, CLASS, SUBCLASS, and RCHI2DIFF fields in the upcom- 
ing Data Release 9. This analysis uses Z.NOQSO redshifts for targets se- 
lected with CLASS_NOQSO="GALAXY" and ZWARNING_NOQSO=0. 
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0.7: at lower redshifts the BOSS LOWZ sample is more dense, 
and we wish to remove overlap between samples. At z > 0.7, the 
efficiency with which we recover redshifts decreases, potentially 
leading to increased systematic errors (Ross et al. 2012). 



3 CATALOG CREATION 

3.1 Target photometry 

Target galaxies are selected as described in Section [2~T] based on 
the best reduced photometry available at the time of target selec- 
tion. All imaging data used by BOSS is based on photometry from 
SDSS Data Release 8 (DR8; [Aihara et al.|20li) . During the early 
phases of BOSS, the final DR8 imaging data were still being pro- 
cessed, and therefore some of the SDSS imaging (9 per cent) used 
for targeting CMASS galaxies is now designated as secondar)ljin 
the DR8 database. Although the measured object parameters from 
different observing runs over the same region agree within the pho- 
tometric errors, there can be significant differences between target 
samples selected from these different observing runs. These dif- 
ferences between our target catalog and that obtained using DR8 
primary photometry arise due to the stochastic variations one ex- 
pects given the magnitude errors and the different photometry. We 
have therefore produced a "combined photometry target" sample 
that uses the photometry input to each run of the targeting soft- 
ware. Thus, rather than thinking of the set of BOSS CMASS galax- 
ies as being a unique sample of galaxies chosen with the properties 
described in Section |2.1| we should really consider the stochastic 
nature caused by photometric errors: we are simply observing one 
of the samples that could have been selected with these properties: 
using different imaging data we would find a different sample. We 
do not expect this issue to have any impact on the analysis or re- 
sults presented here, as it is stochastic in nature. We only include 
this description for completeness and to aid future uses of these 
data. 



3.2 Close-pair corrections 

The protective sheath around each spectroscopic fibre and the fer- 
rule that holds the fiber in the plug plate has a diameter of 62" 
on the focal plane, so no two objects separated by less than this 
can be observed on a single plate. This means that groups and 
clusters of galaxies with members closer than this apparent sepa- 
ration will be systematically under-sampled, strongly affecting the 
measured small-scale clustering signal if uncorrected. The target- 
ing algorithm has been designed to place a fibre on as many objects 
within tight groups as possible (The selection is random with re- 
spect to galaxy properties other than position on the sky). Where 
there are two plates covering a sector we find that approximately 
25 per cent of the pairs separated by < 62" only have one galaxy 
observed; this fraction reduces to < 7 per cent when a sector is 
covered by three or more plates. An algorithm that corrects for 
these effects on small scales is presented and tested in Guo et al. 
pOTT). On large scales, this procedure is equivalent to upweight- 
ing the galaxies nearest to each unobserved galaxy, and we adopt 
this procedure in the analysis presented here. An alternative would 
have been to upweight all galaxies within the sector to compen- 
sate, which would have better shot-noise properties. However, the 
lost galaxies will predominantly be in groups, and thus may have 



i.e. there is an overlapping observation with higher quality photometry 



different clustering properties than the average galaxy. We there- 
fore accept the subsequent slightly increased shot-noise contribu- 
tion. This close-pair correction weight is denoted w cp throughout 
this paper (see also Section [JIT}. For each target we set w cp = 1, 



and add one to this for each CMASS target within 62" that failed 
to get a fiber allocated. This correction affects ~ 5 per cent of all 
the galaxies, with most of these in pairs with w cp = 2. 

3.3 SDSS-II redshifts 

Accurate redshifts for a subsample of the target galaxies were pre- 
viously obtained within the SDSS-II survey ( Abazajian et al. 2009). 
These galaxies were not re-observed by BOSS. We have redshift 
measurements for 100 per cent of the SDSS-II galaxy subsample 
by definition, and although these lie within the BOSS survey re- 
gion, the angular distribution will systematically differ from that of 
the remaining subsample. We do not try to define a survey mask 
that amalgamates the SDSS-II and BOSS observations because the 
SDSS-II redshifts do not even form a random subsample of the 
BOSS targets based on galaxy properties. Instead, we subsample 
the SDSS-II galaxies to match the sector completeness of the galax- 
ies observed within the BOSS program, where sector completeness 
Cboss is defined within each mask sector as 

Nobs + N cp 



Cb 



N t 



targ 



N k 



GO) 



where JV bs, -/V targ , N cp and iVknown are defined in Section [33] 

We also subsample the galaxies so that including the galax- 
ies with previously known redshifts does not change the fraction 
of close pairs observed in any sector. This task is accomplished by 
calculating, for each sector, the fraction of close pairs within the 
sample of targets from which the known galaxies have been re- 
moved. We then subsample close pairs introduced when including 
the known galaxies, randomly removing either galaxy in a pair until 
the fraction of the introduced close pairs matches that of the sample 
without known galaxies. 

Thus we force the known redshifts to have the same statisti- 
cal properties as the BOSS galaxies within each sector, such that 
the angular distribution of the combined sample follows that of the 
BOSS angular mask. 

3.4 Redshift-failure corrections 

We do not achieve stellar classification or a good redshift deter- 
mination for every spectrum taken. The probability of successfully 
obtaining an accurate spectrum is dependent on the fibre used - 
some fibres have degraded transmission, the optical quality (reso- 
lution) is better for spectra in some regions of the CCD than others, 
and the quality of the sky-subtraction is worse where the resolu- 
tion is lower. Bundles of fibres are generally allocated to similar 
regions on the plates, and thus the failure rate is a strong function 
of position in the field-of-view for each observation. This effect is 
shown in Fig. 3 of our companion paper (Ross et al. 2012), where 
specific regions within each field-of-view are demonstrated to have 
worse-than-average failure probabilities. Our overall redshift suc- 
cess rate is 98.2 per cent, so we lack redshifts for a sufficiently 
small subsample that the effect on the measured clustering signal is 
very small. 

We correct for this minor issue by upweighting the nearest 
(based on an angular search) target object for which a galaxy red- 
shift, or stellar classification, has been successfully achieved. The 
redshift distribution of these nearest neighbours matches that of 
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Property 



NGC 



SGC 



total 



^known 
N sts ,r 
iVfaii 
iV cp 

-^missed 

•''used 
N ohs 

JVtarg 

Total area / deg 2 
Effective area / deg 2 



222538 
3766 
7201 
3751 
14116 
4931 



60 792 
1810 
1771 
1122 
3640 
1911 



283 330 
5576 
8972 
4873 
17756 
6842 



207 246 57 037 264 283 
233490 63 685 297 175 
256 303 71046 327 349 



2635 
2584 



709 
690 



3344 
3275 



Table 1. Basic parameters of t he C MASS DR9 sample, when summed over 

We define N x 



all mask sectors (see Section 



2.2 



^-/sector: 



, N x , and 

the meaning of each Nx is given in the text. In our clustering analyses, 
we only consider galaxies with 0.43 < z < 0.7, which is why iV usc( j < 
(iVg a i + -/Vknown)- We split between the Northern Galactic Cap (NGC) 
and Southern Galactic Cap (SGC) regions, for which we calculate separate 
galaxy and random catalogues. The total area is the sum of the areas of all 
mask sectors passing our completeness cut, and the effective area is the sum 
when we multiply the area of each sector by its completeness, Cboss- 



the full sample (see Fig. 4 of |Ross et al.| |2012), suggesting that 
the anisotropic component of the redshift-failure distribution which 
should be the difference between the two does not depend on red- 
shift. Consequently, upweighting the nearest galaxy with a good 
redshift should match the true large-scale density and correct for 
the spatially dependent redshift failure effects. Further details can 
be found in Section 2.3 of |Ross et al.| ( |2012| >. This redshift-failure 
correction weight is denoted w r { throughout this paper (see also 
Section [3~7^ . As in the case of w cp , we define it to be unity for all 
galaxies and then add one if there is a nearby redshift failure. 



3.5 Summary of target objects 

To summarise, the following outcomes are available for BOSS tar- 
gets that are covered by the survey: 

(i) galaxies with redshifts from good BOSS spectra (we denote 
the number in each sector by Ag a i), 

(ii) galaxies with redshifts from SDSS-II spectra (Aknown), 

(iii) spectroscopically confirmed stars (A/star), 

(iv) objects with BOSS spectra from which stellar classification 
or redshift determination failed (Af a ii), 

(v) objects with no spectra, in a close-pair (N cp ), 

(vi) objects with no spectra, or spectra removed following the 
subsampling discussed in Section[3~3) not in a close-pair (iV m i sse d). 



In the following we define the number of target objects per 
sector 

AW = iVstar + iVgal + ATfail + A" cp + ^missed + Aknown, (11) 

and the number of targets observed per sector 

N ohs = AW + Agal + 7Vfa.il . (12) 

The number of good galaxies used in the analysis per sector 
Au S ed, is less than 7Y ga i + ATknown as we only use galaxies with 
0.43 < z < 0.7. Table [T] gives the total split of galaxies in the 
CMASS DR9 target sample into these categories, where we de- 
fine N x — 5"] . N x , and the areas and weighted areas for the 

j. z — /sectors ■*•' & 



CMASS sample in the Northern Galactic Cap (NGC) and Southern 
Galactic Cap (SGC), and combined as derived from the DR9 data. 

Considering the numbers of galaxies in each category per sec- 
tor, we can define a sector completeness as in Eq. [10] and the 
galaxies with previously known redshifts are subsampled to match 
this completeness, as well as the BOSS-only close-pair fraction as 
detailed in Section [33] The distribution of sector completenesses 
across the BOSS footprint is shown in Fig. [T] To remove sectors 
that have only been partially observed, we only retain sectors with 
Cboss greater than 70 per cent. 

We also make a cut on the total redshift failure within each 
sector. First, we define a redshift completeness by 

V - '' (13) 



(-'red — 



N ob 



N Bi 



Then a sector is removed if it has more than 10 BOSS spectra, but 
fewer than 80 per cent of the non-stellar spectra have good redshift 
measurements (i.e. we remove sectors with AT b s — Ac P > 10 and 
Crcd < 0.8). For these sectors we assume that there was a serious 
problem with the observations. Plate 3698 observed MJD 55501 
is responsible for many redshift failures; it was poor data inadver- 
tently included in DR9 with a CMASS failure rate of 23 per cent. 



3.6 Systematic weights 



|RossetaLl ( |2011} have presented a critical examination of the large- 
scale angular clustering of CMASS target galaxies. They demon- 
strated that the density of stars has a significant effect on the ob- 
served density of galaxies, and this can introduce spurious fluctu- 
ations in the galaxy density field. This effect arises because stars 
have a large-scale power signature in their distribution across the 
sky. Additional potential systematics such as Galactic extinction, 
seeing, airmass, and sky background have also been investigated, 
and all have been found to potentially introduce spurious fluctua- 
tions into the galaxy density field, albeit to varying degrees. These 
non-cosmological fluctuations can be corrected for using a weight- 
ing scheme that minimises these fluctuations as a function of a 
given systematic effect (see Fig. 4 of |Ross et al.|201 1) . 

|Ross et aL] ( |2012| investigated systematic effects on the 3D 
clustering of the DR9 CMASS sample. They found that stellar den- 
sity is the primary source of systematic error, and that computing 
a set of weights based on stellar density and ijibi magnitude alone 
has a similar effect as fitting for all five systematic sources simulta- 
neously. Over-fitting these fluctuations can result in removing cos- 
mological power, if the weights remove what are in truth statistical 
fluctuations. Hence the simplicity of correcting for one systematic 
only, with the added dependence on ifn,2, minimises this risk. This 
approach was tested by making use of mock catalogues. We refer 
the reader to Ro ss et aT] < |2012| > for a detailed study of the effect 
of all weighting schemes, and an analysis of each Galactic hemi- 
sphere separately. Note that Ross et al. (2012) explicitly verify that 
the BAO scale is insensitive to these systematic effects. 

The adopted methodology for computing the angular system- 
atic weights used throughout this paper is as follows. The weights 
are defined as 

w B y S (n B , ifib2) — A + Bris (14) 

where n s is the stellar density. A and B are given by 

A = Aoifua + Ai, (15) 

B = Boifon+Bi, (16) 

where the coefficients A = 4.41, Ai = -0.17, Bo = -1.36e-3 
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and B\ — 6.65e — 5 were fitted so as to give a flat relation between 
galaxy density and n s - For ifua < 20.45, A and B were fixed 
at the ifib2 = 20.45 values. These weights were applied to each 
galaxy individually, according to the stellar density of the patch 
of sky in which it lies, and to its observed ifn,2. The stellar density 
map was computed using a HEALPix (Gorski et al. 2005) grid with 
Nside = 256, which splits the sky into equal area pixels of 0.0525 
deg 2 . This pixel size is much smaller than the scale at which the 
systematic effect of stars becomes important (0 > 1°), but large 
enough that the mean number of stars in a pixel is greater than 80 
(implying any shot-noise effects will be small). As this is a large 
scale effect, a relatively coarse mask is sufficient. 



3.7 Final weights and effective volume 

As described in the previous sections, galaxies are weighted to al- 
low for close-pair corrections with w cp , redshift failures with w T { 
and angular systematic weights with w sys . We also apply weights 
to optimise our clustering measurements in the face of shot-noise 
and cosmic variance (Feldman, Kaiser & Peacock 1994), 

Wfkp = 1 , _, - D , (17) 

1 +n(zi)P 

where fi(zi) is the mean density at redshift Zi and Pq — 
20 000/i _3 Mpc 3 . This ignores the scale dependence of the power 
spectrum and chooses a value optimised for the BAO feature, 
Pq ~ P(k = 0.1/iMpc -1 ). We make this simplification for con- 
venience; using the full scale-dependent weights proposed by Feld- 
man, Kaiser & Peacock (1994) does not change our results and 
errors. Ross et al. (2012) find this weighting reduces the variance 
on the CMASS DR9 mock galaxy sample £(s), typically by 20 per 
cent relative to no weighting. Combining wfkp with the systematic 
weights w S ys, the redshift-failure weights w T {, and the close-pair 
weights Wcp, the final weights applied to the galaxies are given by 



lUtot = w FK pW sys (w T { + w cp - 1). 



(18) 



Here wfkp and w sys are mutliplicative weights depending on spa- 
tial location, while w T ! and w cp are additive weights, with default 
of unity. Using these weights, we calculate the effective volume 
using our fiducial cosmology as 



k« = e 



n(zj)P 
l+n(zi)P 



AV(zi 



(19) 



where AV(zi) is the volume of the shell at z; (accounting for the 
observational area). We find V c g = 2.2 Gpc 3 for our CMASS sam- 
ple which covers the redshift range 0.43 < z < 0.7. The weighted 
mean redshift of galaxies within the sample is z e g = 0.57, which 
we define as the effective redshift of our clustering measurements. 



3.8 Random catalog generation 

The evaluation of the correlation function and of the power spec- 
trum requires an estimate of the average galaxy density. To provide 
such an estimate, we generate random catalogs of unclustered ob- 
jects with the detailed redshift and angular selection functions of 
the sample, accounting for the complex survey geometry. In par- 
ticular, the random catalogs account for the differences between 
the Northern and Southern Galactic Caps. To minimise the shot- 
noise induced on clustering measurements, these catalogues have 
70 times more objects than the corresponding galaxy catalogues. 
Numerous tests have confirmed that the survey selection function 



can be factorised into angular and redshift pieces (Ross etal.|2012} . 
The redshift selection function can be taken into account by dis- 
tributing the objects of the random catalogue according to the ob- 
served redshift distribution of the sample. We use the "shuffled" 
catalogue as defined in Ross et aTT| ( |2012| ), where the redshifts are 
matched to randomly selected galaxies. We do this separately for 
the NGC and SGC samples (see Appendix [A] for a further discus- 
sion of this). The completeness on the sky is determined from the 
fraction of target galaxies in a sector for which we obtained a spec- 
trum, with the sectors being areas of the sky covered by a unique 
set of spectroscopic tiles (see Section |2~2} , We upweight close pairs 
and redshift failures in the galaxy catalog as described in Sec- 
tion [34] and therefore include these targets when calculating the 
completeness for the random catalog. Thus the random catalogue 
was subsampled to the sector completeness as given by Eq.|10| 



4 ANALYSIS 

We analyse the BAO feature and fit for a distance to z — 0.57 us- 
ing both the correlation function (Section [51 and power spectrum 
( Section [6l of the 3D galaxy distribution. The steps in both analy- 
ses parallel one another: (i) density-field reconstruction of the BAO 
feature, (ii) computation of the two-point statistics, (iii) estimation 
of errors on these measurements by analysing mock catalogs and 
(iv) extraction of a distance measurement by fitting the data. This 
section details these steps. Details specific to each method as well 
as the results are discussed in later sections. 



4.1 Reconstruction 

As described in Section [T] the statistical sensitivity of the BAO 
measurement is limited by non-linear structure formation. Follow- 
ing Eisenstein et al. (2007a) we apply a procedure to reconstruct 
the linear density field. We emphasise that this improvement is not a 
deconvolution of the correlation function, but uses information en- 
coded in the full density field. In addition to undoing the smoothing 
of the BAO feature, reconstruction also removes the expected bias 
(< 0.5 per cent) in the BAO distance scale that arises from the same 
second-order effects that smooth the BAO feature, which simplifies 
analyses. Reconstruction has recently been applied to the SDSS-II 
DR7 LRG sample atz — 0.35 (Padmanabhan et al. 2012a), and our 
implementation is very similar; we refer the reader there for details 
and simply summarise the steps here: 

(i) Smooth the observed density field to suppress the effects of 
shot-noise and highly non-linear features. We use a Gaussian of 
width I — 15 h~ l Mpc, but demonstrate that our results are insen- 
sitive to this particular choice (see Appendix |Bl[ >. 

(ii) Embed the observed density field into a larger volume with 
a constrained Gaussian realisation. The correlation function of the 
constrained realisation is chosen to match the observed unrecon- 
structed correlations, but we find that our results are insensitive to 
the details of this choice. 

(iii) Estimate the displacements q from the galaxy density field 
<5 ga i using the continuity equation V ■ q = — 5 ga i/& ga i where 6 ga i is 
the galaxy bias. In detail, the above continuity equation is modified 
to account for linear redshift space distortions although we find our 
results are insensitive to the details of this prescription (see Pad- 
manabhan et al. 2012a for the modified continuity equation). The 
galaxy bias 6 ga i is set to a value estimated from the unreconstructed 
correlation function; Appendix |B 1 [ demonstrates that our results are 
insensitive to errors in this choice. 
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(iv) Shift the galaxies by q. Shift the galaxies by an additional 
—fq s s where s is the redshift direction, and / is the logarithmic 
derivative of the linear growth rate with respect to the scale factor. 
This latter shift corrects for linear redshift space distortions. We 
denote this density field by D. 

(v) Generate a sample of points, randomly distributed according 
to the selection function of the survey. Shift these points by q. 
Note that we do not correct the random points for redshift space 
distortions. We denote this density field by S. 

(vi) The reconstructed density field is defined by the difference 
between the density fields defined by D and S. 

Given the large separation between the data in the Northern 
and Southern Galactic Caps, we run reconstruction on these inde- 
pendently. 



4.2 Covariances 

We estimate the sample covariance matrix for the spherically av- 
eraged correlation function and for the spherically averaged power 
spectrum from the distribution of values recovered from 600 galaxy 
mock catalogs. The galaxy mock catalogs are detailed in Manera et 
[aL 1 2012) , and were generated using a method similar to PTHalos 
(Scoccimarro & Sheth 2002), which was calibrated using a suite of 
N-body simulations from LasDamasrj (McBride et al., in prepara- 
tion), and we were able to recover the clustering of halos at ~ 10 
per cent accuracy. A 10 per cent shift in the amplitude of the co- 
variance matrix corresponds to a 3 per cent shift in the error on the 
measured amplitude, and we would expect shifts of a similar order 
in other measured parameters. We consider that knowing our errors 
to this level of accuracy is adequate for our analyses. 

The method can be summarised as follows: we first generate 
600 matter fields at z — 0.55 using second order perturbation the- 
ory (2LPT) for our fiducial cosmology. We choose a fixed redshift 
for simplicity, which assumes that the evolution over the redshift 
range is small. Modelling the evolution in the mocks would re- 
quire a significantly more complex approach, as both the halo mass 
and the method of populating galaxies would both evolve with red- 
shift. We chose to model z = 0.55 since this corresponds to the 
median (unweighted) redshift of the observed galaxy sample, and 
the galaxy mocks were created by matching smaller scale cluster- 
ing without weighting (see Section |3~6Y Since we impose the same 
sampling as the data, the effective weighted redshift of the mocks is 
still z e ff — 0.57. We expect the difference in the covariance matrix 
between z — 0.55 and z = 0.57 to be small and not a significant 
contribution to the errors of our current results. 

We identify dark matter halos using particles from the mass 
field, but we must calibrate the halo mass from these 2LPT ha- 
los (see Manera et al. 2012). We populate the halos with galax- 
ies using a halo occupation distribution (HOD) of the form de- 
scribed in Zheng et al. (2007), with the exact parameter values de- 
termined to reproduce £(r) from the CMASS DR9 data on scales of 
30 - 80 h ~ * Mpc. Although the 2LPT field does not include strong 
non-linear corrections (important on small scales), our method pro- 
duces mocks that match the clustering of CMASS galaxies well 
into the quasi-linear regime when compared with mocks from N- 
body simulations (which do contain the full non-linear evolution). 

The galaxy mock catalogs are initially constructed in 
2400 h' 1 Mpc boxes and then reshaped to fit the DR9 geometry. 



The mock galaxies include redshift-space distortions, follow the 
observed sky completeness, and are down-sampled to correspond to 
the radial number density of the observed data. This was done sepa- 
rately for the NGC and SGC to model the sampling in the data (see 
Appendix [A|. These mock catalogs have also been used in several 
related studies, such as analysing CMASS DR9 systematics (Ross 
et al. 2012), analysis of the clustering of galaxies through simula- 



tions |Nuza et all 



2012), cosmology constraints on redshift space 



Reid et al.|2012 1, and implications of RSD for 



distortions (RSD, 

non-standard cosmologies |Samushia et al. 2012). They were also 
used when comparing evolution between the CMASS and SDSS-II 
Luminous Red Galaxy samples (Tojeiro et al. 2012a), and making 
growth measurements from this comparison (Tojeiro et al. 2012b). 
For this paper, the measurements from the NGC and SGC mocks 
were combined to provide a full measurement for each of 600 mock 
catalogs, just as is done in the observational sample. The galaxy 
mock catalogs and the derived covariance matrices, which we refer 
to as sample covariance matrices, will be made publicly available"] 

The true covariance matrix depends on cosmological parame- 
ters as well as the treatment of galaxy bias, both of which we ne- 
glect. However, we expect this dependence to be relatively weak 
in the parameter range allowed by our data. For example, Labatie 
|et al.| j2012| > study the effect of cosmology dependence in the co- 
variance matrix for the galaxy correlation function for determining 
the acoustic scale; they find that the best-fit value undergoes a small 
shift of 0.3a with a negligible change in the error bar. The effect for 
the CMASS DR9 data should be even smaller as it covers a larger 
volume and has more constraining power equating to less variation 
in parameters than the survey assumed in Labatie et al.]( |2012| l. We 
therefore consider it reasonable to assume no cosmological depen- 
dence in the covariance matrix, which we calculate from galaxy 
mock catalogs based on a fixed, fiducial cosmology. 

Finally, we explicitly test the sample covariance matrix using 
two alternate methods to estimate covariances: 

• a smooth Gaussian model covariance where parameters are fit 
to the galaxy mock catalogs \Xu et al. 2012), and 

• an analytic estimate generated from the monopole power spec- 
trum {deTutter£taQ2012). 

Both of these methods have the advantage of being smooth esti- 
mates of the covariance, which eliminates noise that may compli- 
cate the use of the estimated sample covariance matrix. We find that 
all three methods give consistent results within our quoted errors, 
where we tested by re-fitting the BAO using these three estimates of 
the covariance matrix. Further comparison of methods is provided 
in Manera et al. (2012). For the rest of this paper, we use the sample 
covariance matrices derived directly from the mocks to determine 
errors. 

Although we have 600 independent mocks for each of the 
NGC and SGC samples, these were drawn from the same 600 den- 
sity fields, and so are not fully independent when combined. To 
ensure that this is handled correctly, we form our covariance ma- 
trix for the combined sample by averaging two covariance matrices, 
each calculated from 300 independent joint NGC+SGC samples. 



4.3 Fitting a Distance 

The distance redshift relation for a (thin) redshift slice may be char- 
acterised by the distance to the mean redshift of the sample and 
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its derivative. Getting the former wrong dilates all distances in the 
survey and shifts the BAO feature in the angle averaged correla- 
tion function, but it retains the underlying isotropy of the clustering 
signal. Getting the latter wrong distorts the isotropy of the survey 
and induces higher-order moments in the angle dependence of the 
clustering. Since this work is limited to the angle-averaged cluster- 
ing measurements, we are sensitive primarily to the former. These 
measurements constrain the distance combination 



D v = [cz{l + z) 2 D 2 A H- 1 ] l/i , 



(20) 



where Da is the angular diameter distance and H is the Hubble pa- 
rameter, and we assume that the survey extents are much larger than 
the scales of interest. We also assume that the enhanced clustering 
amplitude along the line-of-sight due to RSDs does not alter the 
relative importance of radial and angular modes when calculating 
spherically averaged statistics (as for the SDSS-II LRGs, [Fercival| 
|et al.|2010) . 

We convert our measurements of the power spectrum and cor- 
relation function into a distance measurement by fitting the acous- 
tic feature to an appropriately dilated template while simultane- 
ously fitting for a "nuisance" broad-band shape. This procedure has 



formed the basis of a number of previous analyses (e.g. Eisenstein 



et al.|2 005 Kazin et al. 2 0"T0l IPercivaTeTa l. 2010, Beutler et al 



2011[|Blake et al.|2011a|bl and[Padmanabhan et al. 2012a} and has 
been explicitly tested to be unbiased even in the presence of in- 
correct assumptions about the underlying cosmology and galaxy 
bias ( |Xu et al.|2012| >. The power spectrum and correlation function 
models are dilated according to P t (k/a) and £t(ar) respectively, 
where P t and £ t are the template functions for a fiducial cosmol- 
ogy. In order to account for the effects of non-linear evolution on 
the BAO feature, the BAO in these templates is artificially smeared 
out according to the prescription given in Eisenstein et al. (2007b). 
Any deviation of the true distance-redshift relation from the fidu- 
cial choice is encoded in a and can be related to the distance to the 
weighted-mean redshift of the sample by 



D v /r 3 - a(D v /r 3 ) s 



(21) 



This distance is relative to the sound horizon r s , for which we 
adopt the definition in Eqs. 4 through 6 of Eisenstein & Hu ( 1998). 
As was discussed in Mehta et al. (2012}, the sound horizon in the 
above is merely a proxy for the distance information encoded in 
the BAO features in the correlation function and power spectrum 
and is therefore insensitive to the choice of definition as long as it 
is consistently used when estimating cosmological parameters. The 
effective redshift of the CMASS DR9 sample is z — 0.57 (see Sec- 
3.6} , and we assume a fiducial cosmology as described in Sec- 
TpThis yields fiducial values D v (z = 0.57) = 2026.49 Mpc, 
= 153.19 Mpc and (£V/r s ) fld = 13.23. 



5 THE CORRELATION FUNCTION 

5.1 Measuring the correlation function 

Fig.[3]plots the CMASS correlation functions before and after re- 
construction, with the best-fit model (see below) overplotted. We 



£(r) = 



(22) 



estimate £(r) using the Landy & Szalay 1 1993 1 estimator 

DD - 2DR + RR 
RR 

where DD, DR and RR are suitably normalised numbers of 
(weighted) data-data, data-random and random-random pairs. For 
the case of the reconstructed correlation function, the DR and RR 



pair counts in the numerator are replaced by DS and SS, where S 
represents the shifted random particles. The errors are estimated by 
applying the same procedure to the mock catalogs and construct- 
ing the sample covariance matrix from the 600 realisations of £(r). 
The average correlation function from the 600 mock catalogs is 
presented in Fig. H] The errors in Fig. [3] are from the diagonal of 
the covariance matrix. We caution the reader that these errors are 
highly covariant, and assessing the significance requires analysing 
the full covariance matrix. 



5.2 Fitting the correlation function 

Our correlation function fits are based on the procedure described 
in |Xu et aT] ( 20 1 2 ) . We give a brief summary of the techniques here. 
Our correlation function model is given by 



A(r) 



where 



k 2 dk 

~2^ 



£m(r) = / '^fPm(k)j (kr)e- 



and 

a / \ a l a 2 



(23) 



(24) 



(25) 



In Eq.[24] the Gaussian term has been introduced to damp the os- 
cillatory transform kernel jo (kr) at high-fc to induce better numer- 
ical convergence. The exact damping scale used in this term is not 
important, and we set a — 1 h~ Mpc, which is significantly be- 
low the scales of interest. The A(r) term is composed of nuisance 
parameters ai.2,3 that help marginalize over the unmodeled broad- 
band signal in the correlation function. Such broadband effects in- 
clude redshift-space distortions, scale-dependent bias and any er- 
rors made in our assumption of the model cosmology. These effects 
may bias our measurement of the acoustic scale if not removed. B 
is a multiplicative constant, allowing for an unknown large-scale 
bias. We use a template P m (k) of the form 



Pm(k) = [PUn(fc) - Pnobao(fc)]« 



-fc 2 £^j/2 



+ Pn, 



>(*). 



(26) 



Eisenstein et al. {2007b 1, Here, Pnn(k) is the linear the- 



as given in 

ory power spectrum and P no bao(fc) is the power spectrum with the 
BAO feature erased. The E„i term is used to damp the acoustic os- 
cillations in the linear theory power spectrum, serving to model the 
effects of non-linear structure growth. We fix S n i = 8 ft -1 Mpc in 
our fits to the pre-reconstruction correlation functions and E n i = 
4 /i _1 Mpc in our fits to the post-reconstruction correlation func- 
tions. We normalise the template to the observed or mock correla- 
tion function being fit at r = 50 h^ 1 Mpc, thereby ensuring that 
B 2 ~ 1. These parameters were tuned on our mock catalogs, and 
we explicitly verify that are results are insensitive to these particu- 
lar choices in Appendix |B2| 

The scale dilation parameter a defined in Eq.pTlcaptures our 
distance constraints; a measures the relative position of the acous- 
tic peak in the data versus the model, thereby characterising any ob- 
served shift. If a > 1, the acoustic peak is shifted towards smaller 
scales, and vice versa for a < 1. 

We obtain the best-fit value of a by computing the x 2 
goodness-of-fit indicator at intervals of Aa = 0.001 in the range 
0.8 < q < 1.2, then identify the value of a that gives the mini- 
mum x 2 and take this as our best-fit value. The x 2 as a function of 
a is given by 



X (a.) — [d — rfb(a)] C [d — A(a)], 



(27) 
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Figure 3. The CMASS correlation function before (left) and after (right) reconstruction (crosses) with the best-fit models overplotted (solid lines). Error bars 
show the square root of the diagonal covariance matrix elements, and data on similar scales are also correlated. The BAO feature is clearly evident, and well 
matched to the best-fit model. The best-fit dilation scale is given in each plot, with the \ 2 statistic giving goodness of fit. 
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Figure 4. Average of the mock correlation functions before and after recon- 
struction showing that the average acoustic peak sharpens up significantly 
after reconstruction. This indicates that, on average, our reconstruction tech- 
nique effectively removes some of the smearing caused by non-linear struc- 
ture growth, affording us the ability to more precisely centroid the acoustic 
peak. 



where d is the measured correlation function and m(a) is the best- 
fit model at each a. C is the sample covariance matrix, and we use 
a fitting range of 28 < r < 200 hT 1 Mpc. We therefore fit over 44 
points using 5 parameters, leaving us with 39 degrees-of-freedom 
(dof). Assuming a multi-variate Gaussian distribution for the fitted 
data (this is tested and shown to be a good approximation in Manera 
|et al.|2012), the probability distribution of a is 



/ \ -x 2 (a)/2 

p(a) oc e k " . 



(28) 



The normalisation constant is determined by ensuring that the dis- 
tribution integrates to 1. In calculating p(a), we also impose a 15 
per cent Gaussian prior on log(a) to suppress values of a < 1 
that correspond to the BAO being shifted to the edge of our fit- 
ting range at large scales. The sample variance is larger at these 



scales, and the fitting algorithm is afforded some flexibility to hide 
the acoustic peak within the larger errors. 

The standard deviation of this probability distribution serves 
as an error estimate on our distance measurement. The standard 
deviation a a for the data and each individual mock catalog can be 
calculated as a\ = (a 2 ) — {a} 2 , where the moments of a are 



(a") = / dap(a)a r ' 



(29) 



Note that (a) refers to the mean of the p(a) distribution in this 
equation only. 

In reference to the mocks, (a) will denote the ensemble mean 
of the a values measured from each individual mock, and a will 
denote the median. The term "Quantiles" will denote the 16* /84 
percentiles, which are approximately the lcr level if the distribution 
is Gaussian. The scatter predicted by these quantiles suffers less 
than the rms from the effects of extreme outliers. 



5.3 Results 

Using the procedure described in |5.2| we measure the shift in the 
acoustic scale from the CMASS DR9 data to be a - 1.016±0.017 
before reconstruction and a = 1.024 ± 0.016 after reconstruction. 
The quoted errors are the a a values measured from the probabil- 
ity distributions, p(a). Plots of the data and corresponding best- 
fit models are shown in Fig. [3] for before (left) and after (right) 
reconstruction. We see that for CMASS DR9, reconstruction has 
not significantly improved our measurement of the acoustic scale. 
However, in the context of the mock catalogues, this result is not 
surprising. 

Fig. [5] shows the a a values measured from the mocks before 
reconstruction versus those measured after reconstruction from the 
correlation function fits. The CMASS DR9 point is overplotted as 
the black star and falls within the locus of mock points. However, 
we see that before reconstruction, our recovered a a for CMASS 
DR9 is much smaller than the mean expected from the mocks. For 
typical cases, reconstruction improves errors on a, but if one has a 
"lucky" realisation that yields a low error to begin with, then recon- 
struction does not produce much improvement. The mock catalog 
comparison in Figure B] shows that the BOSS DR9 data volume 
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Figure 5. Comparisons of <j a errors in mock catalogs before and after re- 
construction as measured from £(r), Reconstruction tends to improve our 
ability to measure a; on a mock-by-mock basis, the average amount of im- 
provement in a a is a factor of 1 .54. However, the amount of improvement 
varies, and 26 (out of 600) of the mocks actually see cr a increase from pre- 
reconstruction to post-reconstruction. The CMASS DR9 point is overplot- 
ted as the black star and falls within the locus of the mock points. 44 (out of 
600) of the mocks have a ratio of cr a after reconstruction compared to be- 
fore reconstruction that is greater than the CMASS DR9 value. Hence, the 
fact that the error on a measured from CMASS DR9 does not decrease sig- 
nificantly after reconstruction is not unexpected in the context of the mocks. 
One can also see that most of the extreme outliers in a a before reconstruc- 
tion have significantly smaller errors after reconstruction. 



is just such a "lucky" realisation, with a strong and well defined 
acoustic peak, and it is therefore unsurprising that reconstruction 
does not reduce the error on a. 

The BAO detection in the CMASS DR9 data is highly sig- 
nificant as illustrated in Fig. |6J Here we have plotted A\ 2 — 
X 2 ~ Xmin> where Xmin 1S the minimum \ 2 that corresponds to the 
best-fit value of a. The dashed line overplotted corresponds to fits 
to the data using a model without a BAO signature. This figure cap- 
tures two tests of BAO significance: the first requires a comparison 
between the solid and dashed curves, and indicates how confident 
we are that the BAO feature exists in the CMASS DR9 data. The 
second uses the plateau height of the Ax 2 curve to indicate how 
confident we are that we have measured an acoustic feature. 

The panel on the left corresponds to our pre-reconstruction 
results and the panel on the right corresponds to our post- 
reconstruction results. Before reconstruction, the minimum of the 
solid curve lies beyond a A\ 2 of 25 from the dashed curve, indi- 
cating that the BAO is detected in CMASS DR9 at greater than 
5rj confidence. Local maxima are seen at greater than A^ 2 of 
36 above the minimum, indicating that the data prefer our best- 
fit value of a at more than 6a. We see similar confidence levels 
post-reconstruction. 

To verify the robustness of our techniques, we also measure 
the best-fit a and a a for our 600 mock catalogues using our fidu- 
cial fitting and reconstruction parameters. We then repeat the same 
fitting with slightly altered models as well as on correlation func- 
tions computed from catalogues that were reconstructed using dif- 
ferent parameters (bias, growth factor and smoothing scale). These 
fitting results are discussed in more detail in Appendix [B] and sum- 
marised in Table IbT] The values in the table are computed after 



discarding the mocks with a a > 0.07. Before reconstruction there 
were 10 such instances, and after reconstruction there was only 1 
such instance. These large uncertainties in the measured a indicate 
that the acoustic signature is weak in these realisations and is there- 
fore not detected with high fidelity (see Xu et al .|2012| for a more 
detailed description of this approach). We find that regardless of 
fitting model parameters or reconstruction parameters, we always 
recover consistent measurements of the a and a a . Hence, our fidu- 
cial model should be trusted to return reliable measurements of the 
acoustic scale. 

Before reconstruction our mocks yield (a) = 1.004 with an 
average error on any single realisation (i.e. the rms or standard de- 
viation) of 0.027 and a standard error on the mean of 0.001. The 
median isa = 1.004 with quantiles of 1 026- After reconstruction, 
we obtain (a) = 1.004 with average error on any single realisation 
of 0.018 and standard error on the mean of 0.001. The median is 
a = 1.004 with quantiles of _o' i8' One can see that given the 
error on the mean, we detect a statistically significant shift in our 
measured mean from the true acoustic scale (a = 1) expected in 
the mocks. This small systematic shift is discussed in more detail 
in Section|7] 

Most important, the average error on a recovered from the 
mocks has decreased after reconstruction. This is illustrated in Fig. 
B] where an overall improvement in a a is evident after reconstruc- 
tion. The greatest improvements occur when the pre-reconstruction 
errors are the worst. The average decrease in a a is a factor of 
1.54, which is equivalent to the effects of increasing the survey 
volume by a factor of 2.3. Therefore, reconstruction appears to sig- 
nificantly improve our ability to measure a precisely, on average. 
This point is further illustrated in Fig.H] where we have plotted the 
average mock correlation function before and after reconstruction. 
One can see the sharpening up of the acoustic peak, indicating the 
effectiveness of the reconstruction algorithm in partially removing 
the smearing of the BAO caused by non-linear structure growth. 
This improvement is what allows a more precise centroiding of 
the peak location. In fitting the average mock correlation function 
before and after reconstruction, we find £„;, the damping of the 
BAO due to non-linear evolution, decreases from 7.58 /iMpc 
to 3.23 /i _1 Mpc. Beyond reducing the distance errors, reconstruc- 
tion also makes our distance estimates more robust to parameter 
choices in our fitting algorithms and reduces the scatter between 
the distance estimates from the the correlation function and the 
power spectrum. We quantify these improvements further in fol- 
lowing sections. 

We next compare the observed scatter in the best-fit a in the 
mocks to the a a estimated in each fit from the X 2 ( a ) curve. In 
Fig. 17] we plot a histogram of (a — {a))/a a from the mocks and 
compare the result to the unit normal distribution. We find excellent 
agreement; a Kolmogorov-Smirnov (K-S) test finds a high likeli- 
hood that the observed distribution is drawn from a unit normal. 
Hence the Gaussian probability distribution obtained from the \ 2 
statistic is an appropriate characterisation of the error on a. 



6 THE POWER SPECTRUM 

6.1 Measuring the power spectrum 

The power spectra recovered from the CMASS DR9 data are shown 
in Fig. [8] before (left) and after (right) reconstruction. The inset 
shows the oscillations in these data, calculated by dividing by a 
smooth model (see Section |6T2| for details). The effect of the re- 
construction algorithm is clear - the large-scale power is decreased 
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Figure 6. Significance of the CMASS DR9 BAO feature before (left) and after (right) reconstruction as measured from f (r). The dashed lines correspond 



to fits to the data using a model without BAO. The quantity A\ 



X 



X-n 



where xini n corresponds to the minimum \ where the best-fit value of 



a lies. By comparing the minimum of the solid curve with the dashed curve, we can quantify how confident we are that the BAO can be measured from the 
CMASS DR9 sample. Similarly, comparing the minimum and the plateau of the solid curve tells us how confident we are that we have measured the correct 
local minima for the acoustic scale. One can see that both before and after reconstruction, we detect the BAO at greater than 5cr confidence and the global 
minimum is itself found within a valley that is 6cr deep. 




for Finger-of-God effects. Given a weight Wi for galaxy i at loca- 
tion Ti, the overdensity field can be written 
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Figure 7. Histogram of (a — (a))/a a measured from £(r) of the post- 
reconstruction mocks, where (a) is the mean. This quantity is a proxy for 
the signal-to-noise ratio of our BAO measurement. We see that this dis- 
tribution is close to Gaussian as indicated by the near-zero K-S D n . The 
corresponding p-value indicates that we are 90 per cent certain our values 
are drawn from a Gaussian distribution, indicating that the values of a a we 
measure from the \ 2 distribution are reasonable descriptors of the error on 
a measured by fitting £(r). 



corresponding to the removal of RSD effects, with the small-scale 
power being further reduced by the reduction in non-linear power. 
These data represent the most accurate measurement of a redshift- 
space galaxy power spectrum ever obtained. 

Power spectra were calculated using the Fourier method first 
developed by Feldman, Kaiser & Peacock ( 1994), as described in 
|Percival et a l. (2007b) and |Reid e t al. (2010JI. We work in redshift- 
space as if observed recession velocities solely arise from the Hub- 
ble expansion. As we focus on measuring angle-averaged Baryon 
Acoustic Oscillations, we do not convert from a galaxy density field 
to a halo density field as in Reid et al. (2010), or apply corrections 



F(r) 



N 



^ w t S D (r t - r) - (w(r)n(r)) 



where TV is a normalisation constant 

( r 1 1/2 

N = < d 3 r{w(r)n(r)f 



(30) 



(31) 



and (w(r)n(r)} is the expected weighted distribution of galaxies 
at location r in the absence of clustering, and n(r) is the galaxy 
density. The quantity Sn is the standard Dirac-5 function. We do 
not apply luminosity-dependent weights (as applied by Percival et 
|al.|2007b|[Reid et al.|2010) , as we are only interested in the BAO, 
and not the overall shape of the power spectrum. 

We chose to model the expected distribution of galaxies us- 
ing a random catalogue with points selected at the mean galaxy 
density (n(r)}, which is then weighted in a similar manner to 
the galaxies. The calculation of this catalogue was described in 
Section |Jl8| The weights in the random catalogue are then renor- 
malised, and compared with the weights applied to the galaxies so 
that J F(r) dr = 0, thereby matching the total weighted number 
density in galaxy and random catalogues. 

Power spectra are calculated using a 2048 grid in a cubic 
box of length 8000 /i -1 Mpc. This zero-pads the galaxies - the 
minimum and maximum galaxy redshifts of the sample correspond 
to distances of 1170 h~ x Mpc and 1780 ft -1 Mpc, so the galaxies 
form an angular sector of a thick shell within this cube. The Nyquist 
frequency for the Fourier transform is approximately 0.8 h Mpc - , 
which is significantly larg er th an the maximum frequency fitted of 
0.3 h Mpc -1 (see Section 6.2 1. The smoothing effect of the cloud- 
in-cell assignment used to locate galaxies on the grid (e.g. Hock- 
ney & Eastwood 1981 . chap. 5) is corrected, and shot-noise is sub- 
tracted following the assumption that galaxies form a Poisson sam- 
pling of the density field (see Feldman, Kaiser & Peacock 1994 for 
details). The power spectrum is then spherically averaged, leaving 



© 201 1 RAS, MNRAS 000,pJJ33 



14 L. Anderson et al. 



Standard 
CMASS DR9 
best-fit model 
X 2 =81.5 / 59 




-1.5 -1 

log 10 k / h Mpc -1 



a, 
S 



a, 







: ■ i ■ ■ ■ 


Reconstructed 


1 




4 + 




• CMASS DR9 
best-fit model 


- 








/=61.1 / 59 


- 


5 

r 




^x 




_ 






*>*•- 




- 


f 


5 
K 

- 2g 
° ? 








- 


5 



|wv* 








0.1 0.2 


0.3 ^ 






k / h Mpc -1 

I 




i V 



-1.5 -1 

log 10 k / h Mpc -1 



-0.5 



Figure 8. The CMASS DR9 power spectra before (left) and after (right) reconstruction with the best-fit models overplotted. The vertical dotted lines show 
the range of scales fitted (0.02 < k < 0.3 h Mpc -1 ), and the inset shows the BAO within this fc-range, determined by dividing both model and data by the 
best-fit model calculated (including window function convolution) with no BAO. Error bars indicate \fCu for the power spectrum and the rms error calculated 
from fitting BAO to the 600 mocks in the inset (see Section|472"]for details). 



an estimate of the "redshift-space" power, binned into bins in k of 
width 0.04/i Mpc" 1 . 

6.2 Fitting the power spectrum 

We fit the observed redshift-space power spectrum, calculated as 
described in Section[6] with a two component model comprising a 
smooth cubic spline multiplied by a model for the BAO, following 
the procedure developed by Percival et al. (2007a c 2010). The 
model power spectrum is given by 



P(k)m = P(fc)smooth X B m (k/a), 



(32) 



where P(fc) S mooth is a smooth model that fits the overall shape 
of the power spectrum, and the BAO model B m (k), calculated for 
our fiducial cosmology, is scaled by the dilation parameter a as 
defined in Eq. [2T1 The calculation of the BAO model is described 
in detail below. This scaling of the acoustic signal is identical to 
that used in the correlation function fits, although the differing non- 
linear prescriptions in (Eqns [23] & [32} means that the non-linear 
BAO damping is treated in a subtly different way. 

Each power spectrum model to be fitted is convolved with the 
survey window function, giving our final model power spectrum to 
be compared with the data. The window function for this convolu- 
tion is the normalised power in a Fourier transform of the weighted 
survey coverage, as defined by the random catalogue, and is calcu- 
lated using the same Fourier procedure described in Section[6](e.g. 
Percival et al. 2007c). This is then fitted to express the window 
function as a matrix relating the model power spectrum evaluated 
at 1000 wavenumbers, k n , equally spaced in < k < 2 h Mpc - , 
to the central wavenumbers of the observed bandpowers k% : 



P(h)m = Y. w ( k *' k n )P(k n ) m - W(h, 0). 



(33) 



The final term W(ki,0) arises because we estimate the average 
galaxy density from the sample, and is related to the integral con- 
straint in the correlation function. In fact this term is smooth (as 



the power of the window function is smooth), and so can be ab- 
sorbed into the smooth component of the fit, and we therefore do 
not explicitly include this term in our fits. 

To model the overall shape of the galaxy clustering power 
spectrum we use a cubic spline ( |Press et al.|1 992), with nine nodes 
fixed empirically at k = 0.001, and 0.02 < k < 0.4 with 
Afc = 0.05, matching that adopted in |Percival et al.| ( |2007c||2010). 
This model was tested in these papers, but we show in Section |B3J 
that it also provides an excellent fit to the overall shape of the DR9 
CMASS mock catalogues, and that there is no evidence for devia- 
tions for the fits to the data. 

To calculate our fiducial BAO model, we start with a linear 



matter power spectrum P(fc)ii n , calculated using CAMB l Lewis et 



|al.|2000[ l, which numerically solves the Boltzman equation describ 
ing the physical processes in the Universe before the baryon-drag 
epoch. We then evolve using the HALOFIT prescription (Smith 
|et al.|20 03), giving an approximation to the evolved power spec- 
trum at the effective redshift of the survey. To extract the BAO, this 
power spectrum is fitted with a model as given by Eq.[32] where we 
adopt a fixed BAO model (Peh) calculated using the Eisenstein & 
|Hu| ( [T9 98) fitting formulae at the same fiducial cosmology. Divid- 
ing P(fc)ii n by the best-fit smooth power spectrum component from 
this fit produces our BAO model, which we denote Pcamb- 

We damp the acoustic oscillations to allow for non-linear ef- 
fects 



tin 



(Be 



l)e 



- k E ni/ 2 



+ 1, 



(34) 



where the damping scale E n ; is a fitted parameter. We assume 
a Gaussian prior on E n ; with width ±2/i -1 Mpc, centred on 
8.24 h~ x Mpc for pre-reconstruction fits and 4.47 h~ x Mpc for 
post-reconstruction fits, matching the average recovered values 
from fits to the 600 mock catalogs with no prior. The exact width of 
the prior is not important, but if we do not include such a prior, then 
the fit can become unstable with respect to local minima at extreme 
values. 
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We fit over scales 0.02/iMpc -1 < k < 0.3/iMpc^ 1 : these 
limits are imposed because the BAO have effectively died out for 
k > 0.3/iMpc" 1 , and scales k < 0.02 /iMpc -1 are sensitive 
to observational systematics ( [Ross et aI7||2012| l. We bin the mea- 
sured power spectrum in k bins of width 0.004 h Mpc -1 , so 70 
data points are included in the fits. The function P(fc) smoot h de- 
pends on 9 free parameters, the amplitudes of the spline nodes. 
Thus, including a and E n ; we fit to 11 parameters in total, and 
the fit has 59 degrees-of-freedom. Goodness of fit between model 
P(k)nt and data is calculated using the \ 2 statistic. We consider 
intervals of Aq = 0.002 in the range 0.7 < a < 1.3 and, for 
each value of a to be tested, we use the Powell routine (Pres s et al] 
1992), starting from a series of widely separated start points, to find 
the spline node values and £„; that result in the minimum value of 

x 2 . 

For each power spectrum fitted, we have estimated the error 
on the best-fit value of a by considering the A% 2 = 1 interval and 
by integrating over the likelihood surface. These measurements are 
found to match extremely well for all of the fits, suggesting that the 
likelihood is well behaved around the minima. We also consider 
the distribution of best-fit a recovered from the mock catalogues, 
as discussed in subsequent sections. For the results presented, in 
ord er to be consistent we adopted the procedure described in Sec- 
tion i 



5.2 



to make measurements from the % surfaces resulting from 



the power spectrum fits. 



6.3 Results 

We have measured the best-fit a and a a using the procedure de- 
scribed in Section [6T2| for power spectra calculated from each of 
600 mock catalogues and from the CMASS DR9 data, either be- 
fore or after applying the reconstruction algorithm. The maximum 
likelihood solution for the dilation parameter from the full CMASS 
DR9 sample is a = 1.022 ± 0.017 before reconstruction and 
a — 1.042 ± 0.016 after reconstruction. Errors were determined 
from the moments of a calculated by integrating the likelihood sur- 
face. The CMASS DR9 power spectra before (left) and after (right) 
reconstruction are plotted in Fig. [8] compared with the correspond- 
ing best-fit models. Before reconstruction, we find Xmin = 81.5, 
while post-reconstruction this reduces to Xmin — 61.1, with 59 
degrees-of-freedom. 

As for the correlation function measurements, we find that re- 
construction does not significantly improve the measurement of the 
dilation parameter from the CMASS DR9 data. This is consistent 
with the results from the mock catalogs: we have measured and fit- 
ted power spectra calculated for all 600 mock catalogs, and Fig. [9] 
presents a scatter plot of the recovered errors before and after re- 
construction. Although reconstruction improves the fit for the ma- 
jority of the mock catalogs, there are a small number for which 
reconstruction increases the recovered error. Also, we see that im- 
provement is more likely where the pre-reconstruction error is high, 
suggesting that variation in the error recovered from different cat- 
alogues is dominated by the "noise" that reconstruction is able to 
remove. In this plot, the star marks the result from the CMASS DR9 
data, showing that the pre-reconstruction error recovered is signifi- 
cantly smaller than the mean expected from the mocks. Given this 
result, we should not be surprised that reconstruction only has a 
small effect on these data. 

The BAO detection from the CMASS DR9 data is highly sig- 
nificant, with a x 2 difference between best-fit models with and 
without the BAO component being approximately 5a before recon- 
struction, dropping slightly to 4.5<r post-reconstruction. The rela- 
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Figure 9. Comparison of mock a a errors before and after reconstruction 
as measured from P(k). This plot is analogous to Figure p\ obtained from 
our fits to £(r). We see the same overall average factor of 1.54 decrease 
in cr a as seen for f (r). This again indicates that in general, reconstruction 
improves our ability to precisely measure a. On a mock-by-mock basis, 
we see a range of results, with 36 (out of 600) mocks actually having cr a 
increase with reconstruction. We once again see that reconstruction does 
not appear to improve the error on our CMASS DR9 measurement of a, 
and that this is not unexpected in the context of the mocks: the CMASS 
DR9 measurement error lies within the locus of values recovered from the 
mocks, with 61 mocks (out of 600) showing a lower fractional improvement 
than that of the CMASS DR9 data. 



tively small difference between significance before and after recon- 
struction matches the difference in a a discussed previously. The \ 2 
surfaces are shown in Fig. |10| before (left) and after (right) recon- 
struction. Low values of a result in the BAO signal being moved 
to large scales where the cosmic variance error increases, which is 
why these models give comparatively good fits. 

From the 600 mocks, pre-reconstruction we recover a mean 
value of the dilation parameter of (a) — 1.004, with average error 
on any single realisation of 0.029 and standard error on the mean 
of 0.002. Post-reconstruction this reduces to (a) = 1.003, with 
average error on any single realisation of 0.019 and standard error 
on the mean of 0.001. There is therefore evidence for a small sys- 
tematic shift between the true value (a = 1) for the mocks and the 
values recovered. A discussion of the systematic errors associated 
with our measurements is provided in Section 17] The decrease in 
mock rms post-reconstruction demonstrates the positive effect that 
reconstruction has on average. Note that when calculating the above 
mean recovered errors we excluded two mocks pre-reconstruction 
for which a a > 0.07, where the BAO feature was not well recov- 
ered. 

The mean values (er a ) match perfectly with the standard de- 
viation of the recovered a values from the mocks, which give 
((a - if) 1 ' 2 = 0.029 and ((a - if) 1 ' 2 = 0.019 pre- and 
post-reconstruction, indicating that the likelihood is extremely well 
behaved. This is not the case if the damping parameter E n ; is fixed 
at an incorrect value in the model to be fitted to the data: insufficient 
damping results in recovered errors that are too small with respect 
to the distribution, while over-damping leads to over-prediction of 
the errors. As in Section [53] we now test the nature of the distri- 
bution of recovered dilation parameters. Fig.fTTJshows a histogram 
of (a — (a))/<7 a compared with a standard Normal distribution. 
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Figure 10. Significance of the CMASS DR9 BAO feature before (left) and after (right) reconstruction as measured from P(k). This figure is analogous to 
Figureplmeasured from our £ (r) analysis. In our P(k) analysis, before reconstruction we detect the BAO in the CMASS DR9 sample at around 5<r confidence, 
similar to our result for £(r), After reconstruction we see a slight drop in the detection level with respect to the pre-reconstruction result. The global maximum 
is found within a valley whose depth is greater than 6cr. 
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Figure 11. Histogram of (a — {a))/a a measured from P(k) of the post- 
reconstruction mocks. This figure is analogous to Fig.lTlobtained from our 
fits to £(r). We again see a near-Gaussian distribution as indicated by the 
small K-S value. This indicates that the a a values we measure from the x 2 
distribution are reasonable estimates of the error on a measured by fitting 
P(k). Note that round-off accounts for the fact that the quoted D n values 
are the same in this Figure and FigurelT] while the p- values differ slightly. 



As is clearly evident, the data are extremely well matched to the 
Gaussian prediction; this is also indicated by the result of a K-S 
test. 

Finally in this section we consider the average BAO signal re- 
covered from the mock catalogues. For each mock, we divide the 
measured power spectrum by the smooth component of the best-fit 
solution convolved with the survey window function. The average 
of these values over all of the mocks is shown in Fig.|12|both before 
and after applying the reconstruction algorithm. The average effect 
of reconstruction is evident on small scales, with the BAO feature 
being enhanced by this algorithm. Fitting to the mocks without as- 
suming a prior on £„; gives average best-fit values of (E„i) = 
8.24 h' 1 Mpc before reconstruction and (£„;) = 4.47 h" 1 Mpc 
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Figure 12. Average BAO signal calculated by dividing the measured power 
from each of the 600 mocks by the best-fit smooth model (solid symbols 
after reconstruction, open symbols before reconstruction). Clearly recon- 
struction enhances the small-scale BAO, where cosmic variance errors are 
significantly reduced. 



following reconstruction, which shows the extent of the improve- 
ment afforded by this technique. Note that these are systematically 
different from the values of £„; recovered from the correlation 
function fits, which results from the way in which the non-linear 
shape was fitted leading to different effective definitions of £„;. In 
the P(k) fits, a multiplicative correction was used, while for £(r), 
an additive correction was adopted: the £(r) fit required less damp- 
ing as he additive correction already acts to damp the importance of 
the BAO component, while the multiplicative correction for P(k) 
afforded by the free shape is itself multiplied by the BAO model, 
and thus more damping is required. 



7 THE DISTANCE TO Z = 0.57 

We now consider how to combine the power spectrum and correla- 
tion function analyses into one estimate for the cosmological dis- 
tance scale. Before reconstruction, we find a = 1.016±0. 017 from 
the correlation function and a = 1.022 ± 0.017 from the power 
spectrum. After reconstruction, we find a = 1.024 ± 0.016 from 
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Figure 13. Comparison of acoustic scale measurements from £(r) and P(k). The left column shows the pre-reconstruction results and the right column 
shows the post-reconstruction results. The top panels show the values of a measured using £(r) versus those measured using P(k); the bottom panels show 
analogous plots for u a . The mock points are shown in grey and the CMASS point is overplotted as the black star. The black cross marks the median values 
of a or cr a along with their quantiles. One can see that there is notable scatter between the values of a and a a measured from the two different statistics. For 
example, a from § and P vary by 0.014 after reconstruction. 
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Figure 14. Comparison of values of a recovered from the mocks before and after reconstruction as measured from £(r) (left) and P(k) (right). As expected, 
there is a tight correlation between the measurements before and after reconstruction, with a slopeshowing the reduced scatter after reconstruction. The 
CMASS DR9 measurements (shown by the stars) lie within the locus of values recovered from the mocks, and the changes in best-fit values seen before and 
after reconstruction are not unusual. 
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Table 2. Key a measurements from BAO in the CMASS DR9 sample. 

a x 2 /dof D v /r s (z = 0.57) 

Before Reconstruction 
f(r) 1.016 ±0.017 30.53/39 13.44 ±0.22 

P(k) 1.022 ±0.017 81.5/59 13.52 ± 0.22 

After Reconstruction 
£(r) 1.024 ±0.016 

P(k) 1.042 ±0.016 

Consensus 1.033 ±0.017 



34.53/39 


13.55 ±0.21 


61.1/59 


13.78 ±0.21 





13.67 ±0.22 
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Figure 15. Histogram of averaged a values from £(r) and P(k) as mea- 
sured from the post-reconstruction mocks. This distribution is very close to 
Gaussian. The near-zero K-S D n value and the corresponding p-value in- 
dicate we are 96 per cent certain that our a values are drawn from a Gaus- 
sian. This result justifies our using a Gaussian probability distribution for 
our CMASS distance measure based on the standard deviation of this dis- 
tribution. 



the correlation function and a — 1.042 ± 0.016 from the power 
spectrum. These measurements are summarised in Table |B2| These 
are small differences, but they approach la and hence demand a 
choice to be made. 

Importantly, our analysis of the mock catalogs shows that this 
level of scatter is not unusual. Fig. [T3] compares the a and errors 
recovered from £(r) and P(k) in our mocks. While the results are 
clearly correlated, there is a notable amount of scatter: about 2. 1 
per cent before reconstruction and 1.4 per cent after reconstruc- 
tion, when one considers the 16-84 per cent quantile. The observed 
small difference in our CMASS measurements before reconstruc- 
tion is unusually good; the larger difference after reconstruction is 
still common, only 1.2<j. Note that here we have not discarded any 
mocks with weak acoustic signals (i.e. a a > 0.07) as we are only 
comparing how £ and P results fare against each other, and not 
examining details of the BAO. Fig.[l4]compares a before and af- 
ter reconstruction for both estimators. Again we find that the shifts 
seem for the CMASS measurements and not unusual given the re- 
sults from the mocks. 

The mocks indicate that both estimators are reasonably unbi- 



ased: the mean a recovered is shifted by only 0.4 per cent from the 
input value, and some of that shift is the actual physical shift of the 
acoustic scale due to non-linear gravity and galaxy clustering bias 
in the mocks (e.g., Crocce & Scoccimarro 2008, Padmanabhan & 



|White|2009HSeo et al.|2010||Mehta et al.|2012) . Furthermore, nei 
ther estimator performs notably better. As shown in the lower pan 
els of Fig 



13 



the errors formed from \ by the two estimators are 
comparable and correlate well on a mock-by-mock basis. 

We therefore conclude that both the correlation function and 
power spectrum estimations of the acoustic scale are appropriate 
and unbiased, but that they include the noise from small scales and 
shot noise differently. Although the power spectrum and correlation 
function are Fourier transform partners in an ideal case, in practice 
this is not true: both functions are considered only over a limited 
domain and with binning. 

Our choice for a consensus distance scale is therefore to aver- 
age the two results. In our cosmological results, we use only the re- 
constructed case. Reconstruction is expected to improve the acous- 
tic scale shifts from non-linear structure formation and galaxy clus- 
tering bias, and it does decrease the scatter in a in our mock cata- 
logs. The actual CMASS data show little change in recovered error 
on a under reconstruction, but this is within the range of behaviours 
of the mocks and is not an argument for avoiding the method. 

To estimate the error bars on the averaged a estimator, we use 
the rms scatter of the results from the mocks for this estimator. Fig. 
T3]shows the distribution of average a values from the mocks. The 
small K-S D n and p-value of 0.96 indicate that we are 96 per cent 
certain our measured a values follow a Gaussian distribution. Since 
we expect our DR9 CMASS a measurement to be Gaussian, using 
the rms of our mock a values as our CMASS error estimate is valid. 
The mocks find a scatter of 1.7 per cent on the average a, which is 
a small decrease from the scatter of 1 .8 per cent on a from £(r) and 
1.9 per cent on a from P(k). This scatter is comparable to the 1.6 
per cent error estimated from X 2 ( a ) °f me CMASS data from the 
fits in both £ and P, which would be another reasonable approach 
to adopt an error. We note that since the averaging does produce 
a small improvement in errors, it is mildly conservative to neglect 
any improvement beyond the errors on the individual estimators. 

We expect the systematic errors in this measurement to be 
much smaller than the statistical errors. The marginalization over 
broadband nuisance terms in the correlation function and over an 
arbitrary broadband spline in the power spectrum gives our tem- 
plate fits great stability to variations in the template or the possibil- 
ity of systematic errors in the galaxy catalog. In Appendix [B] we 
investigate a wide range of variations in the fitting methodologies 
and reconstruction choices, finding variations in a of 0.2 per cent 
or less in all physically reasonable cases. It is more difficult to test 
the effects of systematics in the galaxy catalog, but if we ignore 
all of the corrections for the detected angular systematic variations 
and repeat the clustering analysis, the recovered a value changes 
by only 0.1 per cent despite a notable increase in the large-scale 
power. This result is not surprising: systematic errors of this form 
tend to produce smooth changes in the power spectrum and hence 
are captured by the nuisance terms that we remove in our fitting 
methods. 

|Seo et al.| ( |2008) and |Xu et aT] j2012> show that mild alter- 
ations in the cosmology used for the template in the fit changes 
the recovered Dy jr a at a negligible level, < 0.001 for variations 
consistent with WMAP, when averaged over a number of mock cat- 
alogues. Cosmologies more exotic than the usual cold dark matter 
families of course could open up more dramatic changes; in such 
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cases, one should plan to repeat the fits both to the CMB and BAO 
data sets. 

Our fitting to the mocks does return a value of a that is 0.004 
higher than the input value. As noted above, non-linear structure 
formation and galaxy bias do shift the acoustic scale. |Seo et aL] 
(2010) find shifts of order 0.002 from non-linear structure forma- 
tion, while Meh ta et al.| ( |20lT) finds a similar level from galaxy 
clustering bias. Perturbation theory calculations by Padmanabhan 
& White (2009) yield similar results. However, reconstruction has 
been found to remove these shifts, both in periodic box simulations 
l |Seo et al.|[20T0l |Mehta et al.||20T"fl > and in SDSS-II mock cata- 
logs ([Padmanabhan et al. 2 012a[ ). It is possible that the BOSS DR9 
survey geometry is not large and contiguous enough to remove the 
shifts in full, but it is also possible that the shift in the real data 
might be different than that in the mocks. As the shift is small, we 
have decided not to subtract it from our fitted values and instead to 
consider it as a small systematic uncertainty. 

More exotic galaxy bias models could in principle add addi- 
tional shifts. However, the only physically motivated model known 
that does couple to the acoustic scale is that of Tseliakhovich & 
Hirata (2010), in which relative velocities between the baryons and 
dark matter at high redshift modulate the ability of the smallest ha- 
los to trap gas. Whether this modulation will affect the properties 
of galaxies a million times more massive is speculative. |Yoo et aL] 
pOTT) discuss how the imprint on the acoustic scale in galaxy clus- 
tering could be detected and removed using the three-point func- 
tion, but we have not yet investigated this in the CMASS sample. 

In summary, our consensus value for the acoustic scale fit is 
a — 1.033 ± 0.017. Our estimates of systematic errors are signifi- 
cantly smaller than the statistical error and are negligible in quadra- 
ture. Our best value corresponds to a distance constraint of 



D v (0.57)/r s = 13.67 ±0.22. 



(35) 



We adopt this as our primary result and use it for all of our cos- 
mological interpretations and comparisons to other work. For easy 
reference, the key values of a are summarised in Table [2] For the 
fiducial sound horizon of 153.19 Mpc, Eq. l |35[ ) corresponds to 
L>v(0.57) = 2094 ±34 Mpc. 



8 THE BAO DISTANCE LADDER 

8.1 Comparison to Previous BAO Measurements 

In the last few years, acoustic scale results have been obtained with 
a variety of data sets over a considerable range of redshift. We now 
focus on the comparison between our CMASS DR9 results and past 
work. 

First, we compare the correlation function at z — 0.57 from 
CMASS with that obtained at z = 0.35 by the reconstruction anal- 
ysis of SDSS-II LRGs presented in Padmanabhan et al.| j2012afr. 



Fig. 16 shows these two correlation functions as r £(r). The two 



samples involve different average masses of galaxies and redshifts, 
and hence have a different amplitude of clustering, leading to a 
vertical offset. Both correlation functions use our fiducial ACDM 
cosmology. Given this choice of distance-redshift relation, one can 
see that the acoustic peaks are in excellent agreement. 

Fig.fTTlshows combined significance of the acoustic peak de- 
tection in £(r). In combining the constraints on CMASS DR9 with 
the SDSS-II LRG DR7 data, we neglect the slight overlap in effec- 
tive volumes when using these data in cosmological constraints. 
The LRG data frorrPadmanabhan et al.| (|2012a) cover only the 
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Figure 16. The correlation function measured from CMASS data (black 
circles) versus that from SDSS-II LRG data (grey squares) as shown in 
|Padmanabhan et aT] (2012a) . The vertical offset is due to the difference 
in galaxy bias between the samples; on average the SDSS-II LRGs are 
more luminous and reside in more massive halos. These two analyses used 
slightly different fiducial cosmologies; we have scaled the SDSS-II LRG 
points to cosmology of this paper. One can clearly see that the acoustic peak 
is located at the same position in both datasets. As an aside, we note that 
the difference in the size of the errors has several contributions in addition 
to sample size: the CMASS sample has a higher number density and less 
shot noise, the CMASS sample used 4/i _1 Mpc bins, whereas the SDSS- 
II analysis used 3/i -1 Mpc bins, and the linear scaling of the vertical axis 
causes equal fractional errors to appear larger in the higher bias sample. 




Figure 17. The total significance of the BAO feature, combining the 
CMASS and SDSS-II LRG results, both after reconstruction. This figure 
is analogous to Fig.lSland indicates that in the combined CMASS and LRG 
data sets, we have detected the acoustic peak at greater than 6.5<r, with the 
local minima extending to the ~ Scr level. 



NGC which result in 2496 sq deg of overlapping area with CMASS 
over the redshift interval of 0.43 < z < 0.47. We find this is less 
than 9 per cent of the effective volume of our CMASS sample, and 
less than 5 per cent overlap with the LRG effective volume (frac- 
tionally less since the LRG DR7 data cover a larger area). This is 
consequently a good but not perfect assumption. Combining the 
two correlation functions assuming independence rejects models 
without acoustic oscillations at Ax 2 ~ 45 or 6.7 a. Trying to 
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Figure 18. BAO in the power spectrum measured from the reconstructed 
CMASS data (solid circles with la errors, lower panel) compared with un- 
reconstructed BAO recovered from the SDSS-II LRG data (solid circles 
with lcr eiTors, upper panel). Best-fit models are shown by the solid lines. 
The SDSS-II data are based on the sample and power spectrum calculated in 
|Reidetal.H2010) and analysed by |Percival et al.H2010) ; it has been shifted 
to match the fiducial cosmology assumed in this paper. Clearly the CMASS 
errors are significantly smaller than those of the SDSS-II data, and we also 
benefit from reconstruction, reducing the the BAO damping scale. 
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Figure 20. The BAO distance-redshift relation divided by the best-fit flat, 
ACDM prediction from WMAP (O m = 0.266, h = 0.708; note that 
this is slightly different from the adopted fiducial cosmology of this paper). 
The grey band indicates the 1 a prediction range from WMAP iKomatsu 



etal.|2011). I n addition to the SDSS-II LRG data point from 
et al.H2012a), we also show the result from Percival et al.| 



Padmanabhan 



2010) using a 



combination of SDSS-II DR7 LRG and Main sample galaxies as well as 
2dF Galaxy Redshift Survey data; because of the overlap in samples, we 
use a different symbol. The BAO results agree with the best-fit WMAP 
model at the few percent level. If C m /i 2 were 1 a higher than the best- 
fit WMAP value, then the prediction would be the upper edge of the grey 
region, which matches the BAO data very closely. For example, the dashed 
line is the best-fit CMB+LRG+CMASS flat ACDM model from §|9] which 
clearly is a good fit to all data sets. Also shown are the predicted regions 
from varying the spatial curvature to Q.k = 0.01 (blue band) or varying 
the equation of state to w = —0.7 (red band). 



3000 



,2000 



3 



-.1000 
a 900 
■J 800 
< 700 
% 600 
3 500 
400 



— ' 


1 1 1 1 1 1 


1 ' ' ' 1 ' ' 


— 


- 




B0SS Jl ^^ ll *** , 
^^^WiggleZ 


- 




SDSS-II 






- 


/ 6dFGS 


WMAP ACDM 

1 , , , 1 , , 


- 



0.2 0.4 0.6 
Redshift 



0.8 



Figure 19. A plot of the distance-redshift relation from various BAO mea- 
surements from spectroscopic data sets. We plot Dv{z)/r B times the fidu- 
cial r s to restore a distance. Included here are this CMASS measurement, 
the 6dF Galaxy Survey measurement at z = 0.1 iBeutler et al. 201 1 1, the 
SDSS-II LRG measurement at z = 0.35 ^Padma nabhan et al.|2012a||Xu] 
|et al.|2012||Mehta et al.|2012) , and the WiggleZ measurement at z = 0.6 
(Blake et a ,|201 la). The latter is a combination of 3 partially covariant data 



sets. The grey region is the 1 a prediction from WMAP under the assump- 
tion of a flat Universe with a cosmological constant IKomat su et al.|201 1) . 
The agreement between the various BAO measurements and this prediction 
is excellent. 



place the acoustic peak at other nearby locations and particularly 
at smaller scales is rejected at 8 a. 

Fig.[T8]repeats this comparison with the power spectrum from 
the SDSS-II LRG analysis presented m |Reid et aL"lprJTo| > and[Pef] 
|cival et aL] (f2010). This analysis did not use reconstruction, but one 
can see good agreement in the BAO and significant improvement 
in the error bars with the CMASS sample. 



In Fig. 19 we plot Dv(z) constraints from measurements of 
the BAO from various spectroscopic samples. In addition to the 
SDSS-II LRG value at z - 0.35 ( padmanabhan et al.|2012a) and 
the CMASS consensus result at z = 0.57, we also plot the z = 
0.1 constraint from the 6dF Galaxy Survey (6dFGS) (Beutler et al. 
and a z — 0.6 constraint from the WiggleZ survey (Blake 



2011 



et al.|2011a|. WiggleZ quotes BAO constraints in 3 redshift bins, 



but these separate constraints are weaker and there are significant 
correlations between the redshift bins. We choose here to plot their 
uncorrelated data points for 0.2 < z < 1.0. Each data point here is 
actually a constraint on Dv(z)/r s , and we have multiplied by our 
fiducial r a to get a distance. 

As described further in |Mehta et al.H2012) , the WMAP curve 
on this graph is a prediction, not a fit, assuming a flat ACDM cos- 
mology. For each value of Q m h 2 and flth 2 , one can predict a sound 
horizon, and the angular acoustic scale measured by WMAP plus 
the assumptions about spatial curvature and dark energy equation 
of state then provide a very precise breaking of the degeneracy be- 
tween Q,rn and Ho and hence a unique Dv(z)/r s - Taking the la- 
range of Q m h 2 and Q,i,h 2 produces the grey band in Fig.[l9| There 
is excellent agreement between all four BAO measurements and the 
WMAP ACDM prediction. 
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Following Meht a et al.| ( |2012] >, we divide the distance mea- 
surements by the best-fit WMAP prediction (tt m = 0.266, h — 
0.708) to yield Fig. [20] Focusing first on the data points and the 
grey ACDM region, the data points are consistent with the WMAP 
prediction, but tend to lie closer to the 1 a upward trend in WMAP, 
toward higher Ct m h 2 . In other words, if the WMAP value for Q m h 2 
were 1 a higher, then all of the BAO points would be in superb con- 
sistency with themselves and the CMB under a flat ACDM cosmol- 
ogy; recall that the swath of grey models is a nearly one-parameter 
family so all redshifts move together. We also include here the BAO 
measurement from Percival et al. (20101. Because of the overlap in 
sample with the LRG analysis of Padmanabhan et al. (20 12a I, we 
use a different symbol for this measurement. The correlation func- 
tion analysis from Kazin et al. (2010) gives similar agreement. The 
CMASS BAO value is in perfect agreement with the WiggleZ mea- 
surement ( TBlake et al.|2011a| >. The WiggleZ acoustic scale error is 
3.9 per cent (using their constraint on A(z)), so the CMASS DR9 
error of 1.7 per cent represents a 5-fold improvement in the vari- 
ance. 

Also shown in Fig.[20]is how the WMAP prediction changes 
as one varies the assumptions about dark energy and spatial curva- 
ture. For any specific choice of Qk and w(z), WMAP predicts a 
narrow region set by the range of Q. m h 2 andQbh 2 . Here we present 
the case of Qk — 0.01 with a cosmological constant and a flat 
Universe with w — —0.7; both produce large differences from the 
observations. 

There have also been BAO measurements at 2 ~ 0.55 us- 
ing photometric samples (Padmanabh an et al.|2007||Carnero et al.| 
|2012| |Seo et aL]|2012) . These measure the angular diameter dis- 
tance Da (z) rather than Dy(z). Carnero et al.|{2012[ > measured 
(1 + z)DA/r a — 14.7 ± 1.4 at z = 0.55 using the angular cor- 
relation function (Crocce e t al.|20lT 1 of SDSS DR7 imaging data 
(Abazajian et al. 2009 1. Seo et al. i 2012| l on the other hand mea- 
sured (1 + z)D A /r s = 14.18 ± 0.63 at z = 0.54 using the 
angular powe r spectrum (|Ho et al.||2012| > of the SDSS-III DR8 
imaging data 1 Aihara et alT| 201 1 iP Despite the different estima- 
tors of two-point statistics used, both results consistently show a 
larger distance scale than the prediction of the WMAP best-fit by 
lcr and 1.4cr, respectively. To compare these values with spectro- 
scopic measurements, we correct the difference between Da(z) 
and Dy(z) using the H(z) calculated from the fiducial cosmology, 
while translating the percentage error on Da (z) to be the percent- 
age error on Dy(z) . The deviation from the WMAP prediction will 
be reduced due to using the fiducial H(z) during this transforma- 
tion; the correction yields Dy(z — 0.55)/r s = 13.6 ± 1.3 and 
D v {z — 0.54)/r s = 13.22 ± 0.58, respectively. Extrapolating 
these values from z ~ 0.55 to z = 0.57 assuming the fiducial 
cosmology gives Dv{z — 0.57) jr a = 14.0 ± 1.4 for Carnero et 
Ial1([20l2"|i and D v (z = 0.57)/r s = 13.81 ± 0.61 for |Seo et al. 



(2012). Therefore the photometric BAO measurements show an ex 
cellent agreement with D v (z - 0.57)/r s = 13.67 ± 0.22 from 
the CMASS measurement. It is clear that these photometric BAO 
measurements also fall into the general upward trend relative to the 
WMAP prediction. 

Similarly, there have been spectroscopic BAO measurements 
that attempt to separate the line-of-sight and transverse clustering 
so as to measure H(z) and Da (z) separately 1 Okumura et al. 2008 



7 The DR8 measurement used 10 000 square degrees of the sky that in- 
cludes the coverage of the CMASS DR9 sample. Therefore the overlap in 
volume between the two samples is approximately 30 per cent. 
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Figure 21. The BAO distance-redshift relation divided by the best-fit flat, 
ACDM prediction from WMAP (grey band), overlaid with the distance ra- 
tio measurements from the 3-year Supernova Legacy Survey data iConley 
|et al.|[201 1) (blue band). Here, we have fit the SNe data using Markov 
Chain Monte Carlo in a parameter space of wo, w a , and spatial curva- 
ture. From the chain, we measure the mean and standard deviation on 
Dy (z) / Dy (0.57). We normalise to the CMASS value when plotting on 
this graph; the supernovae are providing only a relative distance measure. 
There is excellent agreement between the supernova distance-redshift rela- 
tion and the measurements from all BAO experiments. Also shown is the 
direct Ho value from Riess et al. 1 201 1 1; there is mild tension between this 
measurement and the BAO and SNe data. 



Chuang & Wang ^On) . These measurements are at lower redshift 
and hence not directly comparable to our CMASS result. However, 
the agreement in the recovered cosmological parameters is good. 

In summary, a precise view of the Hubble diagram from 
baryon acoustic oscillations over the range 0.1 < z < 0.6 is tak- 
ing shape. These measurements appear highly consistent with the 
standard cosmological model. 



8.2 Comparison to Supernova and Direct H Measurements 

We next offer further comparisons to the distance-redshift relation 
from Type la supernovae and direct Ho measurements. Type la 
supernovae can be used to measure relative luminosity distances. 
We use the 3-year Supernova Legacy Survey (SNLS3) results from 
|Conley et al.| pOTTl, including their systematic error treatment. 
Comparing the supernova distance-redshift relation to that of the 
BAO requires some procedure to bin the individual data points or 
to fit a model (see Lampeitl et al.i2010.for a discussion of possible 
procedures). Due to the systematic errors and the absolute distance 
offset, combining the supernova data into redshift bins would nec- 
essarily yield correlated results. Moreover, the supernovae results 
constrain the luminosity distance Dl, not the Dv measured by the 
BAO. Dv (z) requires H(z) as well, which is related to a derivative 
of Dl(z), with a mild dependence on spatial curvature. The need 
for this derivative recommends fitting a model instead of binning. 

We opt to fit a parametric model. We run Markov Chain Monte 
Carlo (described in the next section) for a model space including 
spatial curvature and an equation of state w(z) = wo + w a (l — a). 
We use CMB data from WMAP in addition to the SNLS3 data; 
we opt to include CMB measurements so that the dependence on 
spatial curvature remains mild. CMB data alone would have 3 di- 
mensions of significant degeneracies in this parameter space; the 
supernova data will attempt to break these degeneracies. We then 
use the Markov Chain to infer the constraints on ratios of Dv {z)to 
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Dy(0.57). This tells us what the supernova distance-redshift rela- 
tion predicts for the ratio between two BAO measurements, subject 
to the regularisation of the supernova data implied by the paramet- 
ric cosmological model that we have chosen. In effect, we have fit 
a three-parameter distance-redshift relation to the supernova data 
and then used this to infer Dv(z) from the observed distance mod- 
uli. As a technical note, these results will differ slightly from those 
from MCMC that combine BAO and SNe data, because the BAO 
data will limit the exploration of the distance-redshift degeneracy 
space. 

Comparing z = 0.57 to Z — 0.35, we find that the super- 
novae measure LV(0.35)/ZV(0. 57) = 0.6579±0.0063,a l.Oper 
cent inference. From the CMASS measurement of Dv(0.57)/r 3 = 
13.67 ± 0.22, this predicts D v (0.35)/r s = 8.99 ± 0.14 ± 0.09, 
where the first error arrises from the CMASS error and the second 
error is due to the error in the supernova propagation. This predic- 
tion can be compared to the Padmanabhan et al. (2012a) measure- 
ment from SDSS-II LRG of D v (0.35)/r 3 = 8.88 ± 0.17. Hence, 
the ratio of these two BAO measurements agrees well with the su- 
pernova data. 

Similarly, at z — 0.10, we find that the supernovae measure 
D v (0.10) /D v (0.57) = 0.2018 ± 0.0038, a 1.9 per cent infer- 
ence. The combination with the CMASS data would then predict 
D v (0.10)/r 3 = 2.759 ± 0.044 ± 0.052, following the notation 
from the previous paragraph. This can be compared to the 6dFGS 
measurement of 2.81 ±0. 13, where we have scaled from z — 0.106 
to z — 0.1. Again, the ratio of the BAO measurements agrees well 
with the supernova distance scale. 

We present these results graphically in Fig.[2T] Here, we nor- 
malise the Dy(z) from the Markov Chain at z = 0.57 and con- 
sider the mean and 1 a range explored by the chain. Of course, 
one might have chosen to normalise at another redshift; this ver- 
sion presents how well the CMASS BAO data can be transferred 
to other redshifts. One can see the excellent agreement with all of 
the other BAO results. One also sees that the supernova relative 
distance scale is still more constraining than the BAO relative dis- 
tance scale, by a factor of order 2-3. Of course, the supernovae do 
not provide an absolute distance scale; this plot is indicating only 
their constraint on the slope of the distance-redshift relation. In the 
future, we many wish to combine SNe and BAO distances to fur- 
ther constrain the reciprocity, or distance-duality, relation which is 
a generic prediction of any theory of gravity where photons follow 
null geodesies ( TBassett & Kunz|2004||Lampeitl et al.|2010J . 

Finally, considering the constraint on the Hubble constant, the 
supernovae predict a tight relation between Dv(0.57) and 1/Ho- 
We quote this quantity as #o-Dy(0.57)/(0.57c) = 0.844 with a 
2.3 per cent error. Using this result, the CMASS BAO data with a 
sound horizon given by the fiducial cosmological model predicts 
Ho — 68.9 km/s/Mpc, with 1.7 per cent error from the z — 0.57 
calibration and 2.3 per cent error from the supernova transfer to 
Z — 0. This value is in mild tension with the direct measurement 
of Ho = 73.8 ± 2.4 km/s/Mpc using the NGC 4258 maser and 
HST near-IR observations of Cepheid variable stars (Ries s et al.| 
201 1 ). Fig.pTlplots this measurement, but we remind readers that 
the placement of this point assumes the fiducial value of r s , which 
creates a 1 per cent uncertainty not included in the errors. We will 
quantify this point further in the next section, using a full Markov 
Chain. 
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Figure 22. 68 per cent contours for Ho vs f2 m in the ACDM cosmological 
model. The CMASS DR9 BAO data improve our measurements of Ho and 
Q m , and are consistent with the SDSS-II LRG measurements. The dashed 
grey lines are lines of constant Sl m /i 2 using the WMAP7 values and mod- 
ulated by Iff (fl m h 2 = 0.1334±o:oo6s)- 
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Figure 23. 68 per cent contours for Ho vs Ck in the oCDM cosmologi- 
cal model. The BAO data break the geometrical degeneracy in the CMB, 
and the CMASS DR9 measurements are consistent with the SDSS-II LRG 
measurements. 



9 COSMOLOGICAL PARAMETERS 

To explore the implications of these results for the values of cosmo- 
logical parameters, we consider the standard CDM parametrisation 
of the baryon and matter densities {Q m , ^t}, the primordial power 
spectrum slope n s , the optical depth to reionization r, the Hubble 
constant Ho and the amplitude of matter clustering erg. We also 
examine models with a non-zero curvature Ok as well as models 
where the dark energy differs from a cosmological constant with 
an equation of state parameterised by w(a) = iurj + (1 — a)w a 
( |Chevallier & Polarski|[2f3oT] |Linder|[2053) l, where a is the scale 
factor. 

We follow the methodology in Mehta et al. (2012), using the 
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Table 3. List of the datasets used in the MCMC chains for measuring cosmological parameters. 



Dataset Name Description 



References 



CMB WMAP7 data 

LRG SDSS-II Luminous Red Galaxies 

6dFGRS 6dF Galaxy Redshift Survey sample 

CMASS SDSS-III Data Release 9 Constant Mass Sample 

SN 3-year Supernova Legacy Survey compilation (SNLS3) 

Hq Direct Hubble Constant Measurement 



Komatsuet al. 1 20 Hi 
Padmanabhan et al. 1 20 12a I 



etal. 



Beutl 

This paper 



2011 



Conley etal (201 : 
Riess et al. (20TT ■ 



Table 4. The first two columns show the cosmological model and the data set in each case. The remaining columns show the cosmological parameter values 
estimated from the mean of the posterior distribution. The uncertainties, estimated from the second moments of the distribution are in parentheses and are the 
change to the last significant figures. Eg. 0.268(29) represents 0.268 i 0.029. Empty values correspond to the cases in which the parameter is kept fixed to 
its fiducial value, i.e. Jl^ = 0, miq = — 1> w a = 0. 



Cosmological Model Data Sets 



Qmh 2 



n n 



H 
km/s/Mpc 



Ok 



ACDM 


CMB 


0.1341(56) 


0.268(29) 


71.0(26) 










ACDM 


CMB+CMASS 


0.1392(36) 


0.298(17) 


68.4(13) 










ACDM 


CMB+LRG 


0.1362(33) 


0.280(14) 


69.8(12) 










ACDM 


CMB+LRG+CMASS 


0.1384(31) 


0.293(12) 


68.8(10) 










ACDM 


CMB+LRG+CMASS+6dF 


0.1384(31) 


0.293(12) 


68.7(10) 










ACDM 


CMB+LRG+CMASS+SN 


0.1373(30) 


0.287(11) 


69.2(10) 










ACDM 


CMB+LRG+CMASS+SN+6dF 


0.1373(30) 


0.288(11) 


69.1(10) 










oCDM 


CMB 


0.1344(55) 


0.423(175) 


60.0(123) 


0.039(44) 








oCDM 


CMB+CMASS 


0.1340(53) 


0.299(16) 


67.0(15) 


0.008(5) 








oCDM 


CMB+LRG 


0.1333(53) 


0.278(15) 


69.3(16) 


0.004(5) 








oCDM 


CMB+LRG+CMASS 


0.1336(51) 


0.288(12) 


68.1(11) 


0.006(5) 








oCDM 


CMB+LRG+CMASS+6dF 


0.1336(50) 


0.288(12) 


68.1(11) 


0.006(5) 








oCDM 


CMB+LRG+CMASS+SN 


0.1322(51) 


0.284(12) 


68.3(12) 


0.006(5) 








oCDM 


CMB+LRG+CMASS+SN+6dF 


0.1321(50) 


0.284(12) 


68.2(11) 


0.007(5) 






wCDM 


CMB 


0.1342(58) 


0.263(105) 


75.4(138) 




-1.12(41) 




uiCDM 


CMB+CMASS 


0.1358(59) 


0.323(43) 


65.4(60) 




-0.87(24) 




mCDM 


CMB+LRG 


0.1349(57) 


0.285(25) 


69.0(39) 




-0.97(17) 




mCDM 


CMB+LRG+CMASS 


0.1370(58) 


0.294(27) 


68.6(44) 




-0.99(21) 




uiCDM 


CMB+LRG+CMASS+6dF 


0.1363(51) 


0.298(20) 


67.8(31) 




-0.95(15) 




mCDM 


CMB+LRG+CMASS+SN 


0.1399(37) 


0.280(13) 


70.8(18) 




-1.09(8) 




wCDM 


CMB+LRG+CMASS+SN+6dF 


0.1396(37) 


0.282(13) 


70.4(17) 




-1.08(8) 




oujCDM 


CMB+LRG+CMASS 


0.1345(53) 


0.250(42) 


74.1(70) 


0.008(5) 


-1.31(34) 




omCDM 


CMB+LRG+CMASS+6dF 


0.1334(52) 


0.271(31) 


70.5(43) 


0.007(6) 


-1.14(23) 




otoCDM 


CMB+CMASS+SN 


0.1338(53) 


0.280(17) 


69.2(21) 


0.009(5) 


-1.10(8) 




omCDM 


CMB+LRG+CMASS+SN 


0.1337(53) 


0.275(14) 


69.8(18) 


0.007(5) 


-1.09(8) 




oicCDM 


CMB+LRG+CMASS+SN+6dF 


0.1333(52) 


0.276(13) 


69.6(17) 


0.008(5) 


-1.09(8) 





WoWaCDM 


CMB+LRG+CMASS 


0.1377(58) 


0.282(52) 


70.7(68) 




-1.11(51) 


0.18(122)* 


WOWaCDM 


CMB+LRG+CMASS+6dF 


0.1369(55) 


0.292(41) 


68.9(48) 




-1.02(42) 


0.44(113)* 


uiOi» n CDM 


CMB+CMASS+SN 


0.1389(62) 


0.281(17) 


70.3(23) 




-1.07(16) 


-0.85(96)* 


IDOlBaCDM 


CMB+LRG+CMASS+SN 


0.1392(59) 


0.280(14) 


70.6(19) 




-1.08(15) 


0.10(87) 


WoWaCDM 


CMB+LRG+CMASS+SN+6dF 


0.1385(58) 


0.281(14) 


70.2(17) 




-1.08(15) 


0.08(81) 


ouiou> a CDM 


CMB+LRG+CMASS 


0.1347(54) 


0.263(54) 


72.7(79) 


-0.009(6) 


-1.13(54) 


-0.70(139)* 


OU! 10a CDM 


CMB+LRG+CMASS+6dF 


0.1341(53) 


0.284(40) 


69.2(50) 


-0.009(7) 


-0.93(41) 


-0.93(130)* 


OU>0«JaCDM 


CMB+CMASS+SN 


0.1344(54) 


0.280(17) 


69.5(21) 


-0.012(6) 


-0.91(17) 


-1.31(102)* 


ouio«j a CDM 


CMB+LRG+CMASS+SN 


0.1348(53) 


0.277(14) 


69.8(18) 


-0.012(5) 


-0.89(16) 


-1.44(93)* 


OU>0«JaCDM 


CMB+LRG+CMASS+SN+6dF 


0.1343(52) 


0.278(14) 


69.5(17) 


-0.012(5) 


-0.88(15) 


-1.40(94)* 


OUIfjlOaCDM 


CMB+LRG+CMASS+SN+HO 


0.1364(51) 


0.270(12) 


71.1(15) 


-0.010(5) 


-0.93(16) 


-1.46(95)* 


OUIQlOaCDM 


CMB+LRG+CMASS+SN+H0+6dF 


0.1359(50) 


0.270(12) 


70.8(14) 


-0.010(5) 


-0.93(16) 


-1.39(96)* 



1 The CMB+LRG values are taken from|Mehta et al.|(2012). * Datasets allow chains to explore parameter space outside the prior of —3.0 ^ w a ^ 2.0. 
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Figure 24. 68 per cent contours for Ho vs w in the toCDM cosmological 
model. As in Fig. [23] addition of the BAO data break the degeneracy in 
the CMB data. The differences in the two are due to the different redshift 
dependence of dark energy and curvature. 
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Figure 26. 68 per cent contours for Hq vs w in the tuCDM cosmological 
model comparing different datasets. Contrast this with Fig. |25| the smaller 
redshift lever arm of the BAO data makes them less sensitive to variations 
in the equation of state. 
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Figure 25. 68 per cent contours for Ho vs Six in the oCDM cosmologi- 
cal model comparing different datasets. The SNe data are less effective at 
constraining curvature, given its subdominance at low redshifts. 



CosmoMC {Lewis & Bridle 2002) Markov Chain Monte Carlo sam- 
pler to map the posterior distributions of these parameters. Our 
BAO distance constraints are parameterised as described above as 
a measurement on Dv/r 3 at z = 0.57; we augment these with the 
z — 0.35 measurement from Padmanabhan et al. (2012a) as well 
as the 6dF measurement at z = 0.106 (Beutler e t al.|201l) . These 
measurements have very little overlap in redshift and cover differ- 
ent angular patches, and we treat them independently. We do not 
include the WiggleZ measurements (Blake et al. 201 lab) given 
the significant overlap with the sample presented here. However, 
as discussed in the previous section, the WiggleZ measurements 
agree very well with the distances derived in this work. In addi- 
tion to these BAO measurements, we include observations of the 
temperature and polarization fluctuations in the cosmic microwave 
background (CMB) by the WMAP satellite ( |Komatsu et al.|2011| >, 
as well as measurements of the expansion history by the 3-year 
Supernova Legacy Survey (Conley et al.|201l) and local measure- 



Figure 27. 68 per cent contours for u>o vs £7k in the otoCDM cosmolog- 
ical model for CMB+LRG+CMASS+SN (shaded red), CMB+SN (dashed 
blue), and CMB+LRG+CMASS (dashed black) datasets. Note the relative 
orthogonality of the contours - the BAO data are very effective at constrain- 
ing curvature, while the SNe data constrain the equation of state. Combining 
the two yield tight constraints both orQk and wq . 



ments of the Hubble constant by |Riess et al.| ( |20lT| ). We summarise 
the data sets used in Table|3] 

We summarise our estimated cosmological parameters and 
their uncertainties for different assumptions about the background 
cosmology in Table El The discussion below highlights particular 
cross-sections through this space of models and parameters, fo- 
cusing on comparisons between the LRG and CMASS samples as 
well as comparisons between the cosmological constraints from the 
BAO and supernova data. 

The most restricted model we consider (denoted ACDM) is 
a ACDM cosmology with no spatial curvature; the dark energy is 
assumed to be a cosmological constant with w = —1. As is clear 
from Fig. [22] this model is already highly constrained by the CMB 
through a combination of constraints on the physical matter den- 
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Figure 28. 68 per cent contours for wo vs w a in the uiouiaCDM cosmolog- 
ical model for CMB+LRG+CMASS+SN (shaded red), CMB+SN (dashed 
blue), and CMB+LRG+CMASS (dashed black) datasets. We have used a 
prior on w a as follows: —3.0 ^ w a ^ 2.0. Compare the overlaps in this 
case with Fig. |27| the constraints from the BAO and SNe are less comple- 
mentary. 
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Figure 29. 68 per cent contours for Ho vs Q m (top left), wo vs Qk 
(top right), and w a vs Six (bottom left), and w a vs u>o (bottom right), in 
the owowaCViM cosmological model for CMB+LRG+CMASS+SN (solid 
red) and CMB+LRG+SN (dashed black) datasets. We have used a prior on 
w a as follows: —3.0 ^ w a ^ 2.0. 



sity Q, m h z and the distance to the last scattering surface. However, 
the current WMAP data cannot fully separate fl m and h, leading 
to an uncertainty in both of these measurements along the direction 
of constant fl m h n , where n ~ 3 (Percival et al. 2002). Adding a 
single low redshift distance measurement, from either the LRG or 
CMASS data, significantly reduces this uncertainty. The similar er- 
rors of the two BAO distances lead to similar constraints on Ho' 
±1.2 kms^Mpc -1 for the LRG sample and ±1.3 kms^Mpc -1 
for the CMASS sample (a 1.7 per cent measurement). Combining 
these reduces this error to ±1.0 kms _1 Mpc _1 (a 1.4 per cent mea- 
surement), a reduction by ~ v2 (Table U I. 

Allowing the curvature or ui (for a constant equation of state) 
to vary (denoted oCDM and wCDM respectively) opens up a de- 
generacy in the Qk/wq — Ho plane when only the CMB data are 



considered (Figs,|23"]and|24[l. This degeneracy is broken by the in- 
troduction of a single distance measurement, as one might have 
expected from Fig. [20] The larger degeneracy in fi K — Ho and 
the subsequently tighter constraints from the BAO measurements 
compared to Wo — Ho results from the different redshifts at which 
curvature and dark energy become important. The BAO and CMB 
measurements are connected through the sound horizon and curva- 
ture has the dominant effect on this lever arm. This effect is visually 
apparent in Fig. [20] where the effect of curvature is mostly an offset 
in the distance-redshift relation (over the redshifts for which we are 
sensitive), while changing wo results in a non-trivial change to the 
shape of the distance redshift relation. This result also explains the 
difference in the improvement when the LRG and CMASS samples 
are combined. The two samples do not have a wide enough lever 
arm to improve the constraints on Ho in the wCDM case. By con- 
trast, for the oCDM case, the errors in Ho drop by ~ 25 per cent 
from the LRG only case. 

For both these cosmological models, one can also compare 
the constraints from the SN data with those from the BAO data 
as shown in Figs. [25] and [26] The qualitative difference between 
the SN and BAO distance ladders is that while the SN data are 
a regular distance ladder, building out from low redshift to high 
redshift, the BAO are an "inverse" distance ladder, calibrated at the 
CMB and extending down to low redshift. The SN therefore only 
weakly constrain the curvature (Fig. [25} and are more sensitive to 
Wo, with the reverse being true for BAO. This effect is reflected 
in Figs. [25] and [26] The constraints on curvature are significantly 
improved by the BAO data, and they do not improve significantly 
upon the addition of the SN data. For the i«CDM case, while the 
BAO measurements have lower constraining power, their different 
redshift dependence gives them a different degeneracy direction to 
the SN, resulting in improved constraints. These trends are repeated 
when we consider two parameter models of the expansion history: 
owCDM (Fig.[27j and w w a CDM (Fig.[28}, with the error ellipses 
being more orthogonal when the curvature is allowed to vary. 

None of the individual probes are currently sufficiently sen- 
sitive to constrain the combination of J1k, wo and w a - We there- 
fore combine the SN and BAO data to obtain constraints on these 
models (Fig. |29|>. This cosmological model is also the one recom- 
mended by the Dark Energy Task Force (Albrecht et al. 2006) as 
the baseline to compare different dark energy experiments. They 
recommend using the inverse of the area of the 95 per cent error 
ellipse in the wo — w a plane as a "Figure of Merit" (FoM) for the 
experiment. Our results (CMB+LRG+SN+CMASS) yield a FoM 
of 14.4, compared to a FoM of 1 1 .5 (CMB+LRG+SN) reported by 
Mehta et al. (2012); the improvement driven by the inclusion of the 
higher precision BOSS measurement is clear. 

Finally, as was discussed extensively in |Mehta et al.|p012| ), 
the combination of the SN and BAO distances allows one to transfer 
the CMB distance scale down to the local Universe and constrain 
Ho- Fig. [30] demonstrates that the resulting value of Ho is robust 
to changes in assumed cosmological model. While the difference 
between the inferred value of Ho from these data and the direct 
measurement of |Riess et al.H201~l) is not statistically significant in 
these data (~ la), they may be brought into better agreement by 
adding an additional relativistic energy density component equiva- 
lent to 4.26 ± 0.56 neutrino species, instead of the canonical 3.04 
(see e.g. Mehta et al. 2012 for more details on the mechanism). 
Improvements in both data sets in the future will elucidate if the in- 
troduction of new physics is warranted or if the explanation is more 
mundane. 
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Figure 30. Measured values for H using the CMB+LRG+CMASS+SN 
dataset under various cosmological models. This figure shows that we get 
consistent results for Hq , which is slightly smaller than the direct measure- 
ment by |RiessetaT]p01H . 



10 DISCUSSION 

We have presented the first constraints on cosmology and the dis- 
tance scale from the Data Release 9 CMASS galaxy sample of the 
Baryon Oscillation Spectroscopic Survey. Our results are based on 
accurate 3D positions of 264 283 massive galaxies covering 3275 
sq. deg. with an effective redshift z — 0.57. This is the largest 
sample of the Universe ever surveyed at this high density and the 
derived B AO constraints are the most accurate determination of the 
distance scale within the crucial redshift range where the expansion 
of the Universe begins to accelerate. 

The large survey volume and high sampling density of the 
CMASS galaxies allow the detection of the acoustic oscillations 
predicted by theories of the early Universe at very high significance 
(> 5 a). The acoustic signature is seen both as a clear peak in the 
correlation function and a series of "wiggles" in the power spec- 
trum. The measures are highly consistent, and we use both statistics 
in our final results. We determine the statistical significance of our 
measurements using a large number of mock catalogs based on sec- 
ond order Lagrangian perturbation theory ( Manera et a l.|2012| >, and 
test the covariance matrices so derived with two analytic methods. 
Our analysis of the mock catalogs shows that our measurements 
and their errors are not at all unusual, and would be expected given 
our sampling if the underlying cosmology were of the ACDM fam- 
ily. Applying reconstruction (Eisenst ein et al.|[2007a) to the data 
does not significantly improve our measurement of the acoustic sig- 
nature, which is to be expected based on comparison to mock cata- 
logs since the pre-reconstruction error in the CMASS DR9 data is 
smaller than for a "typical" realisation. 

We obtain a distance measurement from the power spectrum 
and correlation function by fitting the acoustic feature to an ap- 
propriately scaled template, while marginalising over variations in 
the broad-band shape. Our results are very robust to the procedure 
employed to marginalize over broad-band power, and indeed the 
configuration- and Fourier-space fits provide consistent constraints 
even though the template form and procedure are quite different. 
The scale parameter, a, relates Dy /r 3 to the value in a fiducial 
cosmology. Since we use angle-averaged statistics in this work, 



the relevant distance measure is Dv = [cz(l + z) 2 D\/ 'H\ , 
and it is measured relative to the sound horizon, r 3 . Since the cor- 
relation function and power spectrum include noise from small- 
scales and shot-noise differently, we average the two determina- 
tions to obtain our consensus result on the distance to z — 0.57, 
Dv/r s = 13.67 ± 0.22, where we use the scatter in the mock 
catalogs as an estimate of the error (1.7 per cent) on this average. 

|Reid et aT] ( |2012| > and Sanchez et al. (20121 use the correla- 
tion function over a wide range of scales to constrain cosmological 
parameters. We find excellent agreement between their results and 
the pure-BAO measurement described here, despite slightly differ- 
ent choices of binning, fit range, etc. This demonstrates that the dis- 
tance information is dominated by the sharp acoustic feature rather 
than the broad-band power (which we have explicitly marginalised 
over in our analysis). 

The BOSS result can be combined with other BAO measure- 
ments to form an "inverse distance ladder" which tightly constrains 
the expansion rate from z ~ 0.1 to z ~ 0.6. The acoustic signature 
measured in BOSS is in excellent agreement with earlier SDSS re- 
sults ( |Percival et al.|2010[|Padmanabhan et al.|2012a) , and the dis- 
tance to z ~ 0.6 is in almost perfect agreement with that inferred 
by WiggleZ ( TBlake et al.|2011a} . In general the independent BAO 
results are all consistent with the same underlying (flat, ACDM) 
cosmology. Even with only a fraction of the survey completed, the 
BOSS constraint is already the tightest distance constraint in the 
ladder (1.7 per cent), with an errorbar 2.3 times smaller at z ~ 0.6 
than the combined, earlier WiggleZ measurements (Blak e~et al.| 
|2011a} . The BAO distance ladder suggests a slightly larger dis- 
tance scale than the best-fit to the 7-year WMAP data, lying closer 
to the 1 a upper limit in WMAP toward higher Q, m h 2 . With this 
slightly higher value of Q m h 2 , the BAO measurements are in su- 
perb agreement with each other and the CMB within the context 
of a flat ACDM cosmology. While SNe do not provide an absolute 
distance, the relative distance scale inferred from SNLS SNe data 
is in good agreement with that inferred from BAO. 

BOSS continues to amass data, and we expect these con- 
straints to tighten significantly, as data will be collected through 
mid-2014. 



11 ACKNOWLEDGEMENTS 

WJP is grateful for support from the UK Science and Technol- 
ogy Facilities Research Council, the Leverhulme Trust, and the 
European Research Council. AJR is grateful to the UK Science 
and Technology Facilities Council for financial support through 
the grant ST/1001204/1. CGS acknowledges funding from project 
AYA2010-21766-C03-02 of the Spanish Ministry of Science and 
Innovation (MICINN). DJE, XX, and KM were supported by NSF 
grant AST-0707725 and NASA grant NNX07AH11G. SH and 
DJS are supported by the US Department of Energys Office of 
High Energy Physics (DE-AC02-05CH11231). NP and AJC are 
partially supported by NASA grant NNX11AF43G. FP acknowl- 
edges support from the Spanish MICINNs Consolider grant Mul- 
tiDark CSD2009-00064. MAS was supported by NSF grant AST- 
0707266. MECS was supported by the NSF under Award No. AST- 
0901965. 

Funding for SDSS-III has been provided by the Alfred P. 
Sloan Foundation, the Participating Institutions, the National Sci- 
ence Foundation, and the U.S. Department of Energy Office of Sci- 
ence. The SDSS-III web site is http://www.sdss3.org/ 

SDSS-III is managed by the Astrophysical Research Con- 



© 201 1 RAS, MNRAS 000,pp3 



BAO in SDSS-III BOSS DR9 galaxies 27 



sortium for the Participating Institutions of the SDSS-III Collab- 
oration including the University of Arizona, the Brazilian Par- 
ticipation Group, Brookhaven National Laboratory, University of 
Cambridge, Carnegie Mellon University, University of Florida, the 
French Participation Group, the German Participation Group, Har- 
vard University, the Instituto de Astrofisica de Canarias, the Michi- 
gan State/Notre Dame/JINA Participation Group, Johns Hopkins 
University, Lawrence Berkeley National Laboratory, Max Planck 
Institute for Astrophysics, Max Planck Institute for Extraterres- 
trial Physics, New Mexico State University, New York University, 
Ohio State University, Pennsylvania State University, University of 
Portsmouth, Princeton University, the Spanish Participation Group, 
University of Tokyo, University of Utah, Vanderbilt University, 
University of Virginia, University of Washington, and Yale Uni- 
versity. 

We acknowledge the use of the Legacy Archive for Microwave 
Background Data Analysis (LAMBDA). Support for LAMBDA is 
provided by the NASA Office of Space Science. 

The reconstruction computations were supported by facilities 
and staff of the Yale University Faculty of Arts and Sciences High 
Performance Computing Center. The MCMC computations in this 
paper were run on the Odyssey cluster supported by the FAS Sci- 
ence Division Research Computing Group at Harvard University. 

Catalogue construction, power spectrum calculations, and fit- 
ting made use of the facilities and staff of the UK Sciama High 
Performance Computing cluster supported by the ICG, SEPNet and 
the University of Portsmouth, and the COSMOS/Universe super- 
computer, a UK-CCC facility supported by HEFCE and STFC in 
cooperation with CGI/Intel. 



REFERENCES 

Abazajian, K. N., et al, 2009, ApJS, 182, 543 (DR7) 

Ade P., et al., 201 1 , A&A, in press (arXiv:1101.2 022 1 

Aihara, H, Allende Prieto, C, An, D., et al. 201 1, ApjS, 193, 29 

Alcock C, Paczynski B., 1979, Nature, 281, 358. 

Albrecht A., et al., 2006, astro, |arXiv:a stro-ph/0609591 1 

Bassett B., Kunz M„ 2004, PRD, 69, 1012305 

Beutler F, et al., 2011, MNRAS, 416, 3017 

Blake C, Collister A., Bridle S., Lahav O., 2007, MNRAS, 374, 

1527 
Blake C, et al., 201 la, MNRAS, 415, 1707 
Blake C, et al., 201 lb, MNRAS, 415, 2892 
Blanton M., et al., 2003, AJ, 125, 2276 
Cabre, A., & Gaztanaga, E. 201 1, MNRAS, 412, L98 
Cannon, R., et al., 2006, MNRAS, 372, 425 
Camera, A., et al., 2012, MNRAS, 419, 1689 
Chevallier M., Polarski D., 2001, Int. J. Mod. Phys. D 10, 213 
Chuang, C.-H., Wang, Y, 2011, [arXiv: 1102:2251] 
Cole S., et al, 2005, MNRAS, 362, 505 
Colless M., et al., 2003, |astro-ph/0306581| 
Conley A., et al, 2011, ApJS, 192, 1 

Crocce, M., & Scoccimarro, R., 2008, Phys. Rev. D, 77, 023533 
Crocce, M., et al., 201 1, MNRAS, 417, 2577 
Dawson K., et al., 2012, in preparation, 
de Putter R., et al., 2012 [arxiv: 1 1 1.6596] 
Doi M., et al., 2010, AJ, 139, 1628 
Doroshkevich A.G., Zerdovich Ya.B., Sunyaev R.A., 1978 Soviet 

Astronomy, 22, 523 
Drinkwater M„ et al., 2010, MNRAS, 401, 1429 
Eisenstein D.J., Hu W., 1998, ApJ, 496, 605 



Eisenstein D.J., Hu W., Tegmark M., 1998, ApJ, 504, L57 

Eisenstein D.J., et al., 2001, AJ, 122, 2267 

Eisenstein D.J., 2003, in Wide-field multi-object spectroscopy, 

ASP Conference Series, ed. A. Dey. 
Eisenstein D.J., et al., 2005, ApJ, 633, 560 
Eisenstein D. J., Seo H.-J., Sirko E., Spergel D. N., 2007a, ApJ, 

664, 675 
Eisenstein D.J., Seo H.-J., White M„ 2007b, ApJ, 664, 660 
Eisenstein D.J., Bennett C.L., 2008, Physics Today, April, p. 44 
Eisenstein D.J., et al., 201 1, AJ, 142, 72 |arXiv: 110 1.15291 
Feldman H.A., Kaiser N., Peacock J.A., 1994, ApJ, 426, 23 
FukugitaM., et al., 1996, AJ, 111, 1748 
Gaztanaga E., Cabre A., Hui L., 2009, MNRAS, 399, 1663 
Glazebrook K, Blake C, 2005, Astrophys. J., 631, 1 
Gorski, K. M., Hivon, E., Banday, A. J., Wandelt, B. D., Hansen, 

F K, Reinecke, M., & Bartelmann, M. 2005, ApJ, 622, 759 
Gramann, M. 1993, ApJ, 405, 449 
Guo, H., Zehavi, I., & Zheng, Z. 2011, |arXiv:1111.6598| 
Gunn J.E., et al, 1998, AJ, 116, 3040 
Gunn J.E., et al, 2006, AJ, 131, 2332 
Gunn J.E., et al., 2012, in preparation. 
Hamilton A.J.S., Rimes CD., Scoccimarro R., 2006, MNRAS, 

371, 1205 
Ho, S., et al., 2012, ApJ submitted, |arXiv:1 20L2T37) 
Hockney R.W., Eastwood J.W, 1988, Computer simulation using 

particles, New York: McGraw-Hill 
H0g, E., et al., 2000, A&A, 355, 27 
Hu W., White M., 1996, ApJ, 471, 30 
Hiitsi G., 2006, A&A, 449, 891 
Jones D.H., et al., 2009, MNRAS, 399, 683 
Kaiser N., 1987, MNRAS, 227, 1 
Kazin E„ et al., 2010, ApJ, 710, 1444 
Kolb R., et al., 2006, Report of the Dark Energy Task Force 
Komatsu E., et al., 2011, ApJS, 192, 18 

Labatie A., Starck J. L., Lachieze-Rey M., 2012, ApJ, 746, 172 
Lampeitl H„ et al., 2010, MNRAS, 401, 2331 
Landy S.D., Szalay A.S., 1993, ApJ 412, 64 
Larson D., et al., 2011, ApJS, 192, 16 
Lewis A., Challinor A., Lasenby A., 2000, ApJ, 538, 473 
Lewis A., Bridle S., 2002, PhRvD, 66, 103511 
Linder E. V., 2003, PhRvD, 68, 083504 

Lupton R., Gunn J.E., Ivezic Z., Knapp G., Kent S., 2001, 'Astro- 
nomical Data Analysis Software and Systems X", v.238, 269. 
Manera M., et al., 2012, MNRAS submitted 
Maraston C, Strmbck G, Thomas D., Wake D. A., Nichol R. C, 

2009, MNRAS, 394, 107 
Maraston C, et al., 2012, in preparation 
Masters, K. L„ et al., 2011, MNRAS, 418, 1055 
Matsubara T, 2008a, Phys. Rev. D, 78, 083519 
Matsubara T, 2008b, Phys. Rev. D, 77, 063530 
Mehta K. T, Seo H.-J., Eckel J., Eisenstein D. J., Metchnik M., 

Pinto P., Xu X., 201 1, ApJ, 734, 94 
Mehta K, et al., 2012, MNRAS submitted, | arXiv: 1202. 0092 1 
Meiksin A., White M., Peacock J.A., 1999, MNRAS 304, 851 
Miller C.J., Nichol R.C., Batuski D.J., 2001, ApJ, 555, 68 
Noh Y, White M., Padmanabhan N., 2009, Phys. Rev. D, 80, 

123501 
Nusser A., Dekel A., 1992, ApJ, 391, 443 
Nuza S.E., et al, 2012, MNRAS submitted, |arXiv: 1202.60571 
Okumura T, et al., 2008, ApJ, 676, 889 
Padmanabhan N„ et al., 2007, MNRAS, 378, 852 
Padmanabhan N., et al., 2008, ApJ, 674, 1217 



© 201 1 RAS, MNRAS OOO.pJM 



28 L. Anderson et al. 



Padmanabhan N., White M., 2008, Phys. Rev. D, 77, 123540 
Padmanabhan N., White M., 2009, Phys. Rev. D, 80, 063508 
Padmanabhan N., Xu X., Eisenstein D.J., Scalzo R., Cuesta A.J., 

Mehta K.T., Kazin E., 2012a, [arXiv:1202.0090| 
Padmanabhan N., et al, 2012b, in preparation 
Parejko J, et al, 2012, in preparation 
Peebles P.J.E., Yu J.T., 1970, ApJ, 162, 815 
Peebles P. J. E., 1989, ApJ, 344, L53 
Peebles P. J. E., 1990, ApJ, 362, 1 
Percival W.J., et al., 2001, MNRAS, 327, 1297 
Percival W.J., et al., 2002, MNRAS, 337, 1068 
Percival W.J., et al., 2007a, ApJ, 657, 51 
Percival W.J., et al., 2007b, ApJ, 657, 645 
Percival W.J., et al., 2007c, MNRAS, 381, 1053 
Percival W.J., et al., 2010, MNRAS, 401, 2148 
Percival W.J., White M., 2008, MNRAS, 393, 297 
Perlmutter S. et al., 1999, ApJ, 517, 565 
Pier J.R., et al., 2003, AJ, 125, 1559 
Press W.H., Teukolsky S,A., Vetterling W.T., Flannery B.P., 1992, 

Numerical recipes in C. The art of scientific computing, Second 

edition, Cambridge: University Press. 
Reid B.A., et al., 2010, MNRAS, 404, 60 
Reid B.A., et al., 2012, MNRAS submitted 
Riess A. G. et al., 1998, AJ, 1 16, 1009 
Riess A. G., et al., 2011, ApJ, 730, 1 19 
Ross A. J., et al., 2011, MNRAS, 417, 1530 
Ross A. J., et al., 2012, MNRAS submitted 
Samushia L., et al., 2012, in prep. 

Sanchez, A. G., Baugh, C. M., Angulo, R. 2008, 390, 1470 
Sanchez, A. G., et al., 2009, MNRAS, 400, 1643 
Sanchez, A. G., et al., 2012, MNRAS submitted 
Schlafly E.F., Finkbeiner D.P, 201 1, ApJ, 737, 103 
Schlegel D.J., Finkbeiner D.P, Davis M., 1998, ApJ, 500, 525 
Schlegel D., White M., Eisenstein D., 2009, The Astronomy 

and Astrophysics Decadal Survey, Science White Papers #314 

|arXiv:0902.4680| 
Scoccimarro R., Sheth R.K., 2002, MNRAS, 329, 629 
Seo H.-J., Siegel, E. R., Eisenstein, D. J., White, M., 2008, ApJ, 

686, 13 
Seo H.-J., et al., 2010, ApJ, 720, 1650 
Seo H.-J., et al., 2012, ApJ submitted, [arXiv:1201:2172] 
Smith J.A., et al., 2002, AJ, 123, 2121 
Smith R.E., et al., 2003, MNRAS, 341, 1311 
Stoughton C, et al„ 2002, AJ, 123, 485 
Strauss M.A., et al., 2002, AJ, 124, 1810 
Sunyaev R.A., Zel'dovich Ya.B., 1970, ApSS, 7, 3 
Swanson M.E.C., Tegmark M., Hamilton A.J.S., Hill J.C., 2008, 

MNRAS, 387, 1391 
Tegmark M., et al., 2004, ApJ, 606, 702 
Thomas D., et al., 2012, in preparation 

Tojeiro R., et al., 2012a, MNRAS submitted, |arXi v: 1202. 6241] 
Tojeiro R., et al., 2012b, MNRAS submitted 
Tseliakhovich, D., Hirata, C, 2010, Phys. Rev. D, 82, 083520 
Weinberg, D. H., et al, 2012, |arXiv:1201.2434) 
White M., et al., 201 1, ApJ, 728, 126 
Xu, X., et al., 2012, MNRAS submitted, |arXiv: 1202.00911 
Yoo, J., Dalai, N., Seljak, U., 201 1, JCAP, 7, 18 
YorkD.G., et al., 2000, AJ, 120, 1579 
Zel'dovich Ya.B., Kurt V.G., Sunyaev R.A., 1969 Sov. Phys. 

JETP, 28, 146 
Zheng Z., et al., 2005, ApJ, 633, 791 
Zheng Z., Coil A., Zehavi I., 2007, ApJ, 667, 760 



APPENDIX A: COMPARISONS OF THE NGC AND SGC 

In this paper, we analyse the full sample of the CMASS DR9 galaxy 
catalog combining the Northern Galactic Cap (NGC) and Southern 
Galactic Cap (SGC). We justify this choice in this section, since we 
find no significant differences in the clustering beyond acceptable 
statistical fluctuations. 

The DR9 BOSS footprint contains two disjoint regions: a 
2635 deg 2 region in the NGC and a 709 deg 2 region in the SGC. 
The SGC imaging, on average, is at coordinates with larger Galac- 
tic extinction and was taken under conditions with higher airmass 
and sky background, compared to the NGC. However, Ros s et aT] 
(2012) find that none of these factors has a measurable systematic 
effect on the clustering of DR9 CMASS galaxies. Comparing the 
NGC and SGC, |Ross et al.| ( |2012"| found that the projected number 
density is 3.2 per cent larger in the SGC. ISchlafly & Finkbeiner 
(2011) found offsets in photometry between the NGC and SGC 
data (due to a combination of calibration offsets in the SDSS imag- 
ing data and systematic errors in extinction corrections). Applying 
these offsets to our target selection criteria only removes objects in 
the SGC, so we can apply this correction to our catalogue. Doing so 
reduces the difference in projected number density to 0.2 per cent, 
well within the 2 per cent standard deviation found in the mocks. 
The n(z) in the two hemispheres also appear different. Construct- 
ing a covariance matrix based on the mock n(z), Ross et al. 1 2012 1 
found that 6 per cent of consistent samples exhibit larger NGC to 
SGC variations in the n(z). This fraction was increased to 1 1 per 
cent after applying the Schlafly & Finkbeiner (201 1 ) offsets to the 
sample selection. Our approach is to treat the two regions as hav- 
ing separate selection functions, and to optimally combine the pair 
counts from each sample in order to obtain £(s) measurements for 
the full sample. 

|Ross et aL| ( |2012| > found that the clustering in the NGC and 
SGC is generally consistent, to within 2a. Indeed, we find x 2 = 
56.6 (with 44 data points fitted) when comparing the two regions' 
£(s) measurements in the range 28 < s < 200/t -1 Mpc. The 



greatest differences are found close to the BAO scale, and|Ross et 



aT]p012 1 found no treatment of the data (e.g., applying the Schlafly 



& Finkbeiner 201 1 offsets to the sample selection or applying alter- 
native weighting schemes) that could ameliorate this tension. Be- 
fore reconstruction, we find a — 1.000 ± 0.018 in the NGC and 
a = 1.091 ± 0.030 in the SGC; the BAO position differs by 2.6<r. 
After reconstruction, we measure a — 1.012 ± 0.018 in the NGC 
and a = 1.067 ± 0.035 in the SGC. Hence, after reconstruction, 
the BAO positions differ by only 1.4 a. 

In addition, through analysing our mock catalogues, we find 
no difference between the ensemble properties of the NGC and the 
SGC. From the mock catalogues, we recover a mean (a) = 1.005 
with average error on any single realisation of 0.031 and a median 
a — 1.006 with quantiles to'029 m me NGC before reconstruction. 
In the SGC we find a mean (a) = 0.999 with average error in any 
single realisation of 0.047 and a median a = 0.999 with quan- 
tiles i ' 045 . After reconstruction, we recover a mean (a) = 1.004 
with average error on any single realisation of 0.020 and a median 
a = 1.004 with quantiles ±°° 2 ° in the NGC. In the SGC we find 
a mean (a) = 1.006 with average error on any single realisation 
of 0.041 and a median a — 1.006 with quantiles 1q'c>38- T fle mean 
(a) and median a values of the NGC and SGC mocks are consis- 
tent with each other before and after reconstruction. These results 
also suggest that reconstruction not only improves the precision of 
BAO position measurements, but it also makes the likelihood dis- 
tributions more Gaussian. Accepting that a difference greater than 
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Figure Bl. Correlation function of CMASS DR9 galaxies after reconstruc- 
tion, with different curves corresponding to different input parameters of the 
reconstruction code. There is good agreement between the fiducial choice 
of the parameters (black solid line) and other choices. For each case, we re- 
place either the assumed value of the bias b (red and yellow dot-dash lines), 
the growth rate / (blue and green dashed lines), or the smoothing length of 
the density field (cyan solid line). For reference, the dotted line represents 
the CMASS DR9 correlation function before reconstruction. 



1.4 a happens 16 per cent of the time and that |Ross et al.| (2012) 
find no evidence that the SGC £(s) measurements have significant 
systematic uncertainties affecting the BAO scale, we conclude that 
there is no systematic difference between the two hemispheres and 
therefore base our BAO measurements on the combined sample. 



APPENDIX B: ROBUSTNESS TESTS 

Bl Robustness to Reconstruction Parameters 

In this section, we test the sensitivity of reconstruction to our fidu- 
cial values of the bias 6, the growth rate /, and the smoothing length 
I (used to remove shot noise in the computation of the density field). 
In order to do this, we run reconstruction in a sample of 100 PTHa- 
los mocks assuming the standard values for all except one of these 
parameters. In particular, we tested the effect of assuming a bias of 
b = 1.5 and b = 2.2 (i.e. 20 per cent below and above the average 
recovered bias from the mocks, which was adopted as our standard 
value), and also the effect of assuming a growth rate of / = 0.6 and 
/ = 0.9 (again, 20 per cent below and above our standard value). 
We also tested the effect of choosing a more conservative smooth- 
ing length of I — 20/i _1 Mpc. As both the correlation function and 
power spectrum fits give consistent results for our standard recon- 
struction run (see Section 17), for simplicity in this Appendix we 
only report results measuring the acoustic scale from the correla- 
tion function. 

The results are shown at the end of Table |Bl| On average, dif- 
ferent choices for the values of the reconstruction parameters do not 
bias our distance scale measurements by more than (Aa) < 0.3 
per cent with respect to the measurements using the fiducial values. 
However, even these small biases are not significant considering 
their errors. We conclude that our measurements of the distance 
scale are not sensitive to changes in the fiducial values of the re- 
construction parameters over a wide range of values, and in good 
agreement with the fiducial case. 

We also study the dependence on reconstruction parameters of 
the CMASS DR9 measurements of the distance scale from the cor- 
relation function after reconstruction. We select the same cases we 



studied above for the case of the PTHalos mocks, and run recon- 
struction on CMASS DR9 galaxies for each choice of the parame- 
ters. The results are shown in Table lB2l We find that in all cases the 
distance scale measurements are consistent with the results from 
the fiducial case. It is worth noting that the choice of a galaxy bias 
20 per cent smaller (b — 1.5) than the fiducial case drives the mea- 
surement and the errors above the rest of the cases. The reason 
becomes evident in Figure IbT] which shows the correlation func- 
tions after reconstruction for different values of the reconstruction 
parameters. The shape of the correlation function around the BAO 
peak has been distorted by this particular choice of bias, whereas 
all other choices show results more similar to the fiducial case. This 
is an indication that our estimates for the galaxy bias and the cos- 
mological parameters cannot be completely arbitrary if we want 
to reconstruct the density field accurately. However, for reasonable 
values for these parameters, we do not find a large sensitivity of 
our measurements to these parameters, and find our results to be 
consistent with the fiducial case. 

B2 Robustness of Fitting Algorithm for £ (r ) 

We test the robustness of our correlation function fitting model 
by slightly varying the fiducial model parameters and then re- 
performing the fits to see if we recover consistent values of the 
acoustic scale a. These tests are performed on the mocks as well 
as the CMASS DR9 data. Recall that the fiducial model takes on 
the form given in Equations (23) and ( |25) , where we have taken 
E„; = 8fo -1 Mpc before reconstruction and E n ; = 4 ft -1 Mpc 
after reconstruction. In addition, we specify a fiducial fitting range 
of 28 < r < 200 /i _1 Mpc and use the sample covariance matrix. 
Hence, the fiducial model parameters we alter in performing these 
tests are the order of A(r), the value of E n j, the fitting range, and 
the covariance matrix used. In modifying the form of A(r), polyO 
corresponds to A(r) = 0, poly2 corresponds to a 2-parameter 
A(r) = ai/r 2 + a-z/r, and polyA corresponds to a 4-parameter 
A(r) = a-i/r 2 + a 2 /r + a 3 + a A r. 

The fiducial and tweaked model fit results for 600 mocks are 
shown in Table |B1| We remove mock results with poorly mea- 
sured values of a since a BAO feature was not clearly identified 
(<7 Q > 7 per cent). Nearly perfect 1:1 correlations between the 
values of a are measured from the mocks as shown in Fig. |B2| 
The top two panels of Fig. |B2| show the a values measured us- 
ing the fiducial model plotted against the a values measured us- 
ing a smaller fitting range (50 < r < 200 h^ 1 Mpc) both before 
(left) and after (right) reconstruction. The bottom two panels show 
the corresponding plots for a a . Similar 1:1 correlations are seen 
for most of the other "tweaked" models, implying that our fidu- 
cial model returns unbiased measurements of the acoustic scale. 
The only cases that have larger scatter in the correlations are the 
pre-reconstruction polyO and E n ; = 0/t _1 Mpc cases which is 
not surprising. The prior implies that before reconstruction, there 
is non-negligible broadband smooth signal that may bias our mea- 
surement of the acoustic scale and hence a non-zero form for A(r) 
is required to marginalize over this contribution. The latter implies 
that using a BAO model that does not account for the effects of 
non-linear evolution, which are clearly evident before reconstruc- 
tion, will also bias the measurement of a. After reconstruction, the 
scatter in these cases is greatly reduced as reconstruction partially 
undoes large-scale redshift space distortions and non-linear struc- 
ture growth. 

Similar results for the CMASS DR9 data are shown in Ta- 
ble |B2| In general, our choice of model parameters does not af- 
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Table Bl. Fitting results for various models, found by varying the fiducial fitting model and reconstruction parameters on the ensemble distance scale mea- 
surement from the mocks (i.e. the mean (a) and the median a) as well as the difference in the observed distance scale with respect to the fiducial model on 
a mock-by-mock basis (Aa). The results for the fiducial model, and for different broadband A(r) fitting functions (polyO, poly2, polyA), fitting ranges, and 
non-linear damping S n ; of the acoustic scale, are shown for the correlation function before and after reconstruction. We also present the results of fitting with 
a different covariance matrix (ML) derived based on the technique in Xu et al. 1 20121. For our reconstruction tests, we present the effects of changing the 
fiducial galaxy bias by +20 per cent and —20 per cent (b = 1.5 and b = 2.2), the fiducial growth rate by +20 per cent and —20 per cent (/ = 0.6 and 
/ = 0.9), and the smoothing length to 20 h^ 1 Mpc, a more conservative choice than our fiducial smoothing of 15 h~ 1 Mpc. 



Model 



(a) 



a Quantiles (Aa) 



Aa Quantiles (x 2 }/^°/ 



Before Reconstruction 



Fiducial [/] 1.004 0.027 

Fit with polyO. 0.999 0.026 

Fit with poly2. 1.001 0.027 

Fit with poly4. 1.004 0.027 

Fit between 20 < r < 200 h,- 1 Mpc. 1.001 0.028 

Fit between 50 < r < 200 h- 1 Mpc. 1.005 0.027 

Fit with S nJ -+ 0. 1.000 0.030 

Fit with E ni -» S n ; - 2. 1.002 0.028 

Fit with E n ; -+ S n ; + 2. 1.005 0.028 

Fit using ML covariance matrix. 1.003 0.029 



1.004 
1.000 
1.002 
1.004 
1.003 
1.005 
0.999 
1.003 
1.005 
1.005 



+0.026 
-0.026 

+0.024 
-0.024 

+0.025 
-0.025 

+0.025 
-0.025 

+0.025 
-0.028 

+0.025 
-0.026 

+0.028 
-0.026 

+0.026 
-0.025 

+0.026 
-0.027 

+0.023 
-0.028 



-0.005 


0.009 


-0.004 


-0.002 


0.004 


-0.002 


0.000 


0.001 


-0.000 


-0.002 


0.006 


-0.002 


0.001 


0.003 


0.001 


-0.004 


0.015 


-0.005 


-0.001 


0.005 


-0.002 


0.001 


0.005 


0.001 


-0.001 


0.008 


-0.000 



+0.007 
-0.008 

+0.003 
-0.003 

+0.001 
-0.001 

+0.004 
-0.004 

+0.003 
-0.002 

+0.012 
-0.011 

+0.004 
-0.004 

+0.004 
-0.004 

+0.006 
-0.007 



39.60/39 
42.93/42 
41.24/40 
38.27/38 
41.78/41 
34.20/34 
41.60/39 
39.72/39 
40.11/39 
40.07/39 



After Reconstruction 



Fiducial [/] 

Fit with polyO. 

Fit with poly2. 

Fit with polyi. 

Fit between 20 < r < 200 h' 1 Mpc. 

Fit between 50 < r < 200 h' 1 Mpc. 

Fit with S ni -» 0. 

Fit with E nl -> S ni - 2. 

Fit with Y, nl -> £ n , + 2. 

Fit using ML covariance matrix. 



1.004 


0.018 


1.004 


1.002 


0.018 


1.002 


1.004 


0.018 


1.003 


1.004 


0.018 


1.004 


1.010 


0.017 


1.010 


1.004 


0.018 


1.003 


1.003 


0.019 


1.003 


1.003 


0.018 


1.004 


1.006 


0.018 


1.006 


1.004 


0.019 


1.003 


1.004 


0.019 


1.004 


1.003 


0.019 


1.005 


1.003 


0.018 


1.004 


1.004 


0.018 


1.005 


1.006 


0.019 


1.008 



+0.017 
-0.018 

+0.017 
-0.018 

+0.017 
-0.017 

+0.017 
-0.017 

+0.017 
-0.017 

+0.017 
-0.018 

+0.017 
-0.018 

+0.017 
-0.018 

+0.016 
-0.018 

+0.019 

-0.018 



-0.002 
-0.001 
-0.000 
0.006 
-0.001 
-0.001 
-0.001 
0.001 
-0.000 



0.004 
0.001 
0.000 
0.003 
0.002 
0.003 
0.002 
0.003 
0.004 



-0.002 
-0.001 
-0.000 
0.005 
-0.001 
-0.001 
-0.001 
0.001 
-0.000 



+0.003 
-0.004 

+0.001 

-o.ooi 

+0.000 

-o.ooo 

+0.003 
-0.003 

+0.002 
-0.002 

+0.003 
-0.003 

+0.002 
-0.002 

+0.002 
-0.002 

+0.005 
-0.004 



40.95/39 
45.15/42 
42.53/40 
39.94/38 
47.38/41 
34.55/34 
40.87/39 
40.84/39 
41.62/39 
41.02/39 



Fit to recon. with b — > 1.5. 
Fit to recon. with b — > 2.2. 
Fit to recon. with / — > 0.6. 
Fit to recon. with / — > 0.9. 
Fit to recon. with I -> 20 h^ 1 Mpc. 



+0.016 
-0.022 

+0.015 
-0.023 

+0.017 
-0.021 

+0.015 
-0.022 

+0.016 
-0.023 



0.000 
-0.000 
-0.000 
0.001 
0.003 



0.006 
0.006 
0.002 
0.002 
0.007 



0.000 
-0.001 
-0.000 
0.001 
0.002 



+0.006 
-0.006 

+0.006 
-0.005 

+0.001 
-0.002 

+0.001 
-0.002 

+0.005 
-0.005 



42.56/39 
41.01/39 
40.50/39 
41.09/39 
45.04/39 



i 



Aa : 



«[, 



a 



If] 



where i is the model indicated in the first column. 



2 Note that the error on the mean Aa is \/N smaller than the rms from the mocks quoted in the table, where N is the number of mocks. These much smaller 
numbers would indicate that there is a significant detection of the change in the mean as we change fitting model or reconstruction parameters; however, such 

a small change would not be significantly detected in each mock given the dispersion. 



feet the outcome of the fits. A few cases measure slightly larger or 
smaller values of a, but all fall well within the lcr errorbars. 

We also investigate our measurements of BAO significance 
with respect to the form of A(r). The results are shown in Fig. B3 
after reconstruction. The right panel shows the difference in \ 
between a fit to the data using a model containing BAO and a 
fit to the data using a model without BAO. These curves demon- 
strate how well we have detected the BAO in the CMASS DR9 



data. The solid black curves correspond to subtracting the solid 
line from the dashed line in Fig. [6] The other lines correspond to 
various other forms of A(r), some with more and some with less 
nuisance parameters. Here, the more negative A\ 2 is, the more a 
model containing BAO is preferred. Allowing more or less flexibil- 
ity in the broadband marginalization as parameterised by A(r) does 
not change the fact that a model containing BAO is favoured and 
we have a robust detection of the BAO in the CMASS DR9 data. 
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Figure B3. Ax 2 for CMASS DR9 using various forms of A(r). These plots are analogous to Fig.pl except we have split the two tests of BAO significance 
into separate panels. The left panel shows how robustly we have detected the BAO in the CMASS DR9 sample and the right panel shows how confident we are 
that we have measured the correct acoustic scale. In the left panel, we have plotted the difference in x 2 between 2 fits to the data, one using a model containing 
BAO and one using a model without BAO. We see that this Ax 2 is consistently around —30 for all forms of A(r) indicating that the amount of flexibility 
in the broadband marginalization (i.e. the number of nuisance parameters in A(r)) does not have a significant impact on how well we detect the BAO in the 
CMASS DR9 sample. In the right panel, we have plotted the Ax 2 of the minimum as a function of a. The various forms of A(r) all identify the same best-fit 
value of a and this best-fit is at a Ax 2 well below the plateau in the curve. However, it appears that lower orders of A(r) allow more confident measures of 
a, possibly due to the increased flexibility in higher order forms to fit noise. Regardless, we have at least a 6<r measurement of best-fit a in all cases which is 
robust. 
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Figure B2. Comparison o f res ults obta ined by fitting the mocks using the 
fiducial model, Equations J23) and J25) with a 28 < r < 200 /i _1 Mpc 
fitting range, and those obtained using a 50 < r < 200 h Mpc fitting 
range. The pre-reconstruction results are shown in the left column and the 
post-reconstruction results are shown in the right column. The top panels 
show comparisons of a values and the bottom panels show comparisons of 
a a values. The black cross marks the median values of a or u a along with 
their quantiles. Nearly perfect 1:1 correlations exist between the a values 
measured using the fiducial fitting range and those measured using a slightly 
smaller fitting range both before and after reconstruction. The agreement 
between cr a values is already good before reconstruction, however after re- 
construction, the scatter is reduced to nearly zero. The results of robustness 
tests against other fit parameters are listed in Table |B1| These imply that 
our fiducial model is robust against small changes in model parameters and 
hence should return unbiased measurements of a. 



The actual confidence level changes slightly between the different 
A(r) forms; however, the variation is small and consistently falls 
between 5 — 6a. 

The right panel shows the A\ 2 values from the minimum (or 
best-fit value) and demonstrates how well we have measured the 
acoustic scale. The solid black curve is identical to the solid line in 
Fig.p] The other curves correspond to various other forms of A(r). 
In all cases, the minima lie at the same value of a with the plateaus 
lying at significant Ax 2 above the minima. Although Ax 2 shows 
significant variation between the A(r) forms, we see at least a 6er 
(Ax 2 ~ 36) preference for the best-fit value of a. It appears that a 
lower order or less flexible form for A(r) may return a at a higher 
confidence, which indicates that higher order A(r) may afford the 
model with enough flexibility to start fitting noise. 

We have tested the robustness of the fitting methods by mak- 
ing Gaussian realisations of the correlation function from the full 
covariance matrix and a template for the correlation function. We 
have considered the model defined by Eq. [23] and computed the 
best fit parameters for every simulation. The sample of recovered 
parameter values that we recovered has a Gaussian distribution, 
as expected. We also checked that our simulations are truly Gaus- 
sian by computing the % 2 estimator for each simulation at the "true 
model" (with a = 1, B — 1, and A(r) = 0), and verifying that it 
follows a x 2 -distribution with v = 50 degrees of freedom, which is 
the number of bins in r used in this test. We have found that, while 
in the NGC the inferred errors from the fits agree very well with 
the width of the distribution of a, in the SGC the measured errors 
tend to be slightly under-estimated (by about 0.25 <j a ). 

B3 Robustness of Fitting Algorithm for P(k) 

We have tested that our model for the power spectrum, calculated as 
described in Section [672] provides an adequate match to the power 
spectra of the mocks. In fact the data plotted in Figs.[8l&|12|already 
show this result to some extent as we plot the deviation between the 
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Table B2. Fitting results for various models. Here we explore the effects of 
varying the fiducial fitting model and reconstruction parameters on our mea- 
surements of the distance scale from CMASS DR9. The results for the fidu- 
cial model, for different broadband A(r) fitting functions (polyO, poly2, 
polyA), fitting ranges, and non-linear damping £ n ; of the acoustic scale 
are shown for the correlation function before and after reconstruction. We 
also present the results of fitting with a different covariance matrix (ML) 
derived based on the technique in Xu et al. (2012). For our reconstruction 
tests, we present the effects of changing the fiducial galaxy bias by +20 per 
cent and —20 per cent (b = 1.5 and b = 2.2), the fiducial growth rate by 
+20 per cent and —20 per cent (/ = 0.6 and / = 0.9), and the smooth- 
ing length to 20 h~ 1 Mpc, which is a more conservative choice than our 
fiducial smoothing of 15 h~ 1 Mpc. 



Model 






a 


x 2 


Before Reconstruction 


Fiducial [/] 






1.016 + 0.017 


30.53/39 


Fit with polyO. 






1.018 + 0.020 


40.84/42 


Fit with poly2. 






1.017 + 0.016 


30.74/40 


Fit with polyA. 






1.016 + 0.017 


30.33/38 


Fit between 20 < r 


< 200 h~ 


_1 Mpc. 


1.020 + 0.017 


32.47/41 


Fit between 50 < r 


< 200 h~ 


4 Mpc. 


1.018 + 0.018 


22.99/34 


Fit with £„; -> 0. 






1.005 + 0.013 


30.84/39 


Fit with £ n! -> £„ 


1-2. 




1.012 + 0.015 


29.93/39 


Fit with £„; -> £„ 


1+2. 




1.019 + 0.019 


32.02/39 


Fit using ML covariance matrix. 


1.022 + 0.018 


30.64/39 


After Reconstruction 


Fiducial [/] 






1.024 + 0.016 


34.53/39 


Fit with polyO. 






1.026 + 0.017 


41.82/42 


Fit with poly2. 






1.025 + 0.015 


36.12/40 


Fit with polyA. 






1.024 + 0.017 


33.29/38 


Fit between 20 < r 


< 200 h~ 


^Mpc. 


1.031 + 0.018 


47.31/41 


Fit between 50 < r 


< 200 h- 


_1 Mpc. 


1.022 + 0.016 


25.94/34 


Fit with £ n ; -> 0. 






1.019 + 0.015 


34.18/39 


Fit with £ n! -+ £„ 


l~2. 




1.020 + 0.015 


34.27/39 


Fit with £„; -> £„ 


1+2. 




1.029 + 0.017 


35.10/39 


Fit using ML covariance matrix. 


1.022 + 0.017 


34.30/39 


Fit to recon. with b 


-> 1.5. 




1.033 ± 0.020 


42.97/39 


Fit to recon. with b 


-* 2.2. 




1.021 + 0.015 


46.89/39 


Fit to recon. with / 


-> 0.6. 




1.024 + 0.015 


33.19/39 


Fit to recon. with / 


-> 0.9. 




1.025 + 0.017 


36.53/39 


Fit to recon. with I - 


-> 20 ft" 1 


Mpc. 


1.026 + 0.015 


43.79/39 




Figure B4. Average residual recovered from fitting to the power spectra de- 
rived from the 600 mock catalogues after reconstruction (solid circles with 
lcr errors). The shaded region shows the expected error for any single fit, 
while the solid line shows our fiducial BAO model. Clearly there is no ev- 
idence for any large deviations between the model and data, which might 
have indicated that the spline was unable to match the input power spec- 
trum. 




minimum x 

Figure B5. Histograms of recovered minimum \ 2 values from fits to power 
spectra measured from the mock catalogs: pre-reconstruction (upper panel) 
and post-reconstruction (lower panel). The smooth solid line gives the ex- 
pected distribution of \ 2 values for 59 degrees of freedom, and the vertical 
lines show the values recovered from fits to the CMASS DR9 data. 



measured power spectra and the smooth model: The consistency 
between the data plotted and the expected BAO model shows that 
any residual differences between data and model are of significantly 
lower order than the BAO signal. 

To test the goodness-of-fit further, Fig. |B4| displays the average 
residual recovered after fitting to the 600 power spectra derived 
from the mock catalogs 

{p(h) - p fit (fc0) = ^ E h fc ») - pflt ( fc *)] . < B1 ) 



where P(ki) are the measured band powers, and P (hi) is the 
best-fit model as defined in Eq.|33] As can be seen in Fig. |B4| the 
average residual is well below the scales of both the BAO and the 
difference expected for any single fit (shown by the shaded region 
in Fig. |B4} , so there is no evidence of a systematic inability to fit 
the shape of t he power spectrum over the fitted fc-range 
Fig 



B5 



shows histograms of the recovered best-fit \ values 
from the fits to the 600 mock catalogs, before (upper panel) and af- 
ter (lower panel) reconstruction. These values match the expected 
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distribution of x values for a fit with 59 degrees-of-freedom, 
which is also shown in this plot. This agreement gives us confidence 
that the fit is behaving as expected for the power spectra derived 
from the mock catalogues. If the model was unable to adequately 
fit the mocks, we should expect the recovered x 2 minima to be 
significantly offset from the expected distribution, i.e. if the model 
failed to adequately fit the shape of the recovered power spectrum, 
then we would find systematically worse x 2 values compared with 
those expected. The distribution actually agrees remarkably well 
with that expected, which gives us confidence that the model de- 
scribed in Section [6~2| is adequate for these data. 

The x 2 values from the fits to the data fall within the distribu- 
tion of values from the mocks although, for the pre-reconstruction 
measurement, only 18/600 mocks give a worse x 2 value. However, 
we know from the analysis presented in Section [63] that the pre- 
reconstruction catalogue also gives a smaller-than-average error, so 
this result is perhaps not that surprising. In conclusion we find no 
evidence that the fitting method applied to the power spectra is not 
adequate for recovering the BAO scale. 
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